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GAS RELEASE AND CONDUCTIVITY 


MODIFICATION STUDIES 

Lewis M. Linson and David C. Baxter 

Science Applications, Inc. 

La Jolla, California 92038 

ABSTRACT 

The behavior of gas clouds produced by releases 
from orbital velocity in either a point release or venting 
mode is described by the modification of snowplow equations 
valid in an intermediate altitude regime. Quantitative 
estimates are produced for the time dependence of the 
radius of the cloud, the average internal energy, the 
translational velocity, and the distance traveled. The 
dependence of these quantities on the assumed density 
profile, the internal energy of the gas, and the ratio 
of specific heats is examined. The new feature is the 
inclusion of the effect of the large orbital velocity. 

The resulting gas cloud models are used to calculate 
the characteristics of the field-line-integrated Pedersen 
conductivity enhancements that would be produced by the 
release of barium thermite at orbital velocity in either 
the point release or venting modes as a function of 
release altitude and chemical payload weight. 
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I . INTRODUCTION 

1 2 

A number of experiments ’ have been suggested for 
the Spacelab that require the release of large amounts of 
gas into the ionosphere. Examples of such experiments 
are the: 

• Generation of acoustic gravity waves from a known 
source of momentum and energy; 

• Enhancement of the ionospheric electron concentration 
and modification of the ionospheric conductivity 

by the deployment of easily ionizable vapors such 
as barium or cesium; 

e Dispersement of chemically reactive gases, such 
as H 2 or HgO, in order to affect the neutral and 
ionization chemistry and create holes in the F- 
region of the ionosphere; 

e Creation of large-scale neutral winds by transferring 
momentum to the ambient atmosphere; 

e Testing of Alfvln's critical velocity hypothesis 
regarding the anomalous ionization of neutrals 
moving at high velocity across a magnetized plasma. 

***"•' mm 
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These experiments and others like them require the 
release of large amounts of gas in order to either disperse 
material, transfer momentum, or input energy. The experi- 
ments can be divided into two classes depending on the 
suggested method of deployment — either as a point release 
or in a venting mod 0 We define a point release as a rapid 
or explosive release of gas or material that takes place 
in a time less than a second. In this case, the region of 
space into which the material is dumped is confined and 
the expansion of the gas can be treated as being spherical. 
Other experiments require the deployment of gas in a venting 
mode lasting from five to ten or more seconds. The released 
gas is then spread out in a trail and it expands cylindrically . 

If conducted from Spacelab the released gas initially 
will have a large translational velocity equivalent to a 
hypersonic Uach 10 flow. The large kinetic energy associated 
with the initial translational velocity is a new feature 
that has not been treated previously. This feature 
dominates the interaction with the ambient atmosphere 
leading to rather striking results that affect the properties 
of the release. In order to carry out the objectives of 
the experiments, one needs to be able to describe various 
aspects of the released gas such as what happens to the 
momentum and energy, how fast the cloud expands, how hot 
the expanding cloud becomes, and how far it travels. 



This report considers three aspects of these 

experiments in detail. First, we review the properties 

of the gas cloud expansion that results from a point 

3 4 

release that was obtained from a previous study. ’ 

Second, we develop a model that describes the expansion 
of gas released in the venting mode. Third, we apply 
the results of the above two studies in order to calculate 
the enhancement in conductivity that can be produced by 
the release of barium vapor in sunlight as a function of 
altitude and the amount of release vapor. The major 
part of this report is concerned with the technical 
details of these calculations. 

We have developed models of the dynamic response 
of gas clouds, released in both the point release and 
venting modes, when they interact with the ambient 
rarefied atmosphere at the hypersonic orbital velocity. 
These models are based on a number of simplifying 
assumptions first introduced by Stuart® in order to 
describe the evolution of a stationary point explosion 
in the upper atmosphere. The principal assumptions 
relate to the form of the velocity distribution within 
the cloud and involve a self-similar density distribution 
during the cloud expansion. We have added the effect of 
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the translational velocity at release to Stuart's model 
and have treated the different geometries and time- 
dependencies of the released gas for the two release 
modes. 

The models apply to an intermediate altitude 
regime. This regime is identified by comparing an 
important characteristic cloud size with the molecular 
mean-free-p th, A. For a spherical point release, 
this size is called the equal-mass-radius. a o> and 
is defined as the radius of a spherical volume of ambient 
atmosphere containing a mass equal to the mass of 
released gas. M o< When a Q is comparable to X, neither 
continuum dynamics nor kinetic theory is appropriate for 
describing the interaction of the gas cloud. At lower 
altitudes, where X is much smaller than a o> continuum 
dynamics or fluid mechanics might be applicable. Flow 
of the ambient atmosphere around the gas cloud and the 
effects of shocks and waves created in the ambient 
atmosphere by the gas cloud would be important. On 
the other hand, at higher altitudes where the mean-free 
path is much larger than a Q , a kinetic theory approach 
is required. In this case the majority of the ambient 
molecules would pass right through the roughly spherical 
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gas cloud which is initially traveling at the orbital 
velocity while only a fraction of them would strike one 
of the gas cloud molecules, scattering it out of the 
cloud, thus producing a broad diffuse wake. The 
distinguishing feature of our model is the assumption 
that in the intermediate altitude regime the ambient 
atmosphere interacts strongly with the gas cloud and is 
carried with it, becoming a part of the cloud itself. 

g 

Bernhardt recently completed a similar study 
using a kinetic approach appropriate to the high-altitude 
regime. He used the same criteria for defining the 
high-altitude regime, a Q < 3\, as we adopt in this report. 
Because the physical assumptions appropriate to the two 
regimes are different, the model calculations produce 
qualitatively different results. As one example of this 
difference, the distance that the bulk of an initially 
spherically gas cloud travels is independent of the 
amount of gas released in the high-altitude regime (it 
scales only with X) but in the intermediate-altitude 
regime this distance scales with a Q and thus depends 
on the mass of released gas. 

The parameters that characters© the properties 
of the gas cloud in which we are interested are: the 

7 


< 



scale size, R, of the cloud; the expansion velocity, R, 
of the outer edge of the cloud; the translational velocity, 
V, of the released gas; an equivalent "temperature," T, 
determined by the average internal energy per particle, 
and the distance, L, that the released gas has traveled 
\. 4 .th respect to the ambient atmosphere. In the point 
release case, these parameters are given as functions of 
time since release. The principal difference between 
the point-release case and the venting-mode case is that 
in the frame of the venting canister there is a steady- 
state cylindrically-symmetric distribution of gas down- 
stream of the canister. 

In Sec. II we discuss several topics ” elated to the 
snowplow model as applied to point releases. In Sec. II-A 
we describe the basic assumptions and equations of the 
snowplow model. In Sec. II-B we briefly review and 
summarize the results for the point release of a spherical 
gas cloud at orbital velocity which is treated comprehen- 
sively in Ref. 3. We shall make use of the cloud size 
as a function of time, R(t), in other sections of this 
report and we shall compare the qualitative behavior of 
the gas cloud for this case with similar behavior for 
the ventlng-mode case to be developed in more detail in 
later sections. 
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The neutral gas cloud and wake expansion models 
are based on an assumption that the density distribution 
remains self-similar, i.e., that the cloud maintains 
its density distribution daring the expansion. One 
result of our models is that different density distri- 
butions do not result in qualitative changes in the 
nature of the solutions. In Sec. II-Cwe modify the 

5 

Stuart model to include the fact that at late times 
the gas cloud which initially follows a Stuart snowplow 
expansion eventually evolves into a Gaussian shape 
expanding by diffusion. W>.i treat the effect of diffusion 
as if it were superimposed on the expansion process 
modeled by the self-similar snowplow expansion. The 
result of this analysis is that a cloud with an initial 
constant density profile rapidly evolves into a Gaussian- 
like shape and that the Gaussian scale length of the 
Gaussian distribution has a time dependence almost 
identical to the cloud radius described from the snowplow 
model . 

In Sec. Ill we describe the venting-mode model and 
the equations that determine the cloud parameters. This 
quantitative model conserves the mass, momentum, and 
energy of the expanding and moving gas cloud wake consisting 
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of the released gas and swept -up ambient atmosphere. 

The assumptions on which the model is based and the 

limitations of the model are discussed. Our model 

applies only *o the downstream portion of the wake 

and does not provide a -ascription valid in the 

transition zone between the spherical distribution near 

the venting canister and the cylindrical distribution 

farther downstream. The appropirate conservation 

equations are derived. The characteristics of the vented 

gas cloud rre determined by the model and a specification 

* 

of the parameters associated with the released gas; M Q . 
the rate at which mass exits the canister in kg/s; E 0> 
the rate at which energy exits the canister in J/s; and 
Y, the ratio of specific heats. 

In Sec. IV the results of tho venting model are 
described. The wake-model equations are integrated 
numerically and an accurate analytic approximation is 
derived. In particular, the radius of the wake, the 
radial expansion velocity, the axial velocity of the 
gas in the wake, the temperature, and the distance that 
the gas has traveled with respect to the atmosphere are 
obtained. The dependence of the solutions on the 
various input parameters is illustrated in a series of 


10 



figures. It is found that the analytic expression 
provides a good description of the scaling of these 
solutions with .he input parameters. The wake-model, 
and thus the solutions, becomes more valid l'ar downstream 
at large radius in the wake; the model is less valid at 
small R. For any particular set of the input parameters, 
the wake-model allows- a family of solutions corresponding 
to different assumed initial values of the radial expansion 
velocity, specified at a small value of radius, R Q . 

The approacn of the various numerical solutions to the 
analytic curve at large R is also illustrated by appropriate 
sets of figures. As an example of a quantitative result, 
it is found that far downstream all solutions have the 
behavior R(x) «* Kx* , where the constant K depends on the 
values of the input parameters. 

In Sec. V we apply the results of our theoretical 
modeling of the expansion of the neutral gas clouds to 
the calculation of the change in Pedersen conductivity 
that results then the neutral gas being released is 
neutral barium atom vapor. Much of the motivation for 
the calcMlutions reported here is based on previous 
modeling of large barium icn clouds released in the 
ionosphere from rockets. In particular, the approach 
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used here relies heavily on detailed data analyses of 
the Spruce ion cloud, a barium release experiment 
conducted on 1 February 1971 at 2352 UT at an altitude 
of about 190 km from Kglin AFB in Florida. The data 
consist of both optical and radar data obtained with 
an incoherent scatter radar, an HF long-pulse radar, 
and by radio-wave propagation through the ion clouds. 
These data provide information on the electron concen- 
tration and distribution and scale size of the '.on 
cloud. We combine our recent work on the expansion 
of neutral clouds released from orbital velocity in 
both the point release and venting modes with detailed 
ion cloud modeling performed on the Spruce ion cloud 
in order to calculate the height-integrated Pedersen 
conductivity tnat can be produced by the release of 
various amounts of vapor at a range of altitudes. 

In Sec. V-A we consider the deposition of barium 
ions from the neutral barium atom vapor. Based on 
'he results of Sec. III-C, we assume that the density 
distribution in the barium cloud during the time of ion 
deposition is essentially a Gaussian distribution. We 
make extensive use of earlier work on the deposition of 
iORlzati«.c from barium releases conducted from rockets 
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in order to determine parameters of the barium ion cloud 
in terms of parameters associated with the neutral cloud. 

The procedure that we apply takes into account the competing 
processes of barium oxidation by ambient molecular oxygen 
and photoionization by sunlight. These effects result 
in an altitude-r ^pendent inventory of available barium 
ions from a given amount of barium vapor and in an 
effective ionization time constant associated with the 
loss of neutral barium atoms. The model described here 
also includes the effect of the time delay before 
photoionization begins when the initial density is high 
and impedes the transport of the metastable-exciting 
radiation to the center of the cloud. The choice of 
par&meterizations is based on the detailed analysis 
of the Spruce ion cloud. 

The important parameters for carrying out 
conductivity modification experiments are the field-line- 
integrated Pedersen conductivity and the transverse 
scale size of the ion cloud perpendicular to the 
magnetic field. The scaling of these two parameters 
with altitude and amount of gas released is assumed 
to vary as the radius of the neutral cloud at an appro- 
priate time. The results of the analysis are shown in 
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Sec. V-B in three pairs of figures, each pair corresponding 
to one of the following three cases: a point release 

from a rocket, from Spacelab, and a vented release from 
Spacelab. The figures show contours of constant perpen- 
dicular scale size of the ion cloud and constant Pedersen 
conductivity in the altitude versus chemical-payload- 
weight plane. As a typical example, a 48-kg barium release 
from a rocket at 190 km altitude will produce a field-line- 
integrated Pedersen conductivity of approximately 30 mho, 
several times the ambient ionospheric conductivity, and 
will have a transverse scale size just less than 3 km. 

If the same chemical payload is released at orbital 
velocity, the resulting conductivity is lower by a factor 
of 5 from that produced by a release from a rocket and 
the transverse size of the ion cloud is increased by a 
factor of approximately 2.3. 
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II. REVIEW OF SNOWPLOW MODELS 


The expansion of neutral gas clouds from point 
releasee in the intermediate altitude regime can be 
described by snowplow equations. Point releases at 
orbital velocity are described by a model that represents 
an extension of Stuart's^ original model. The venting 
model described in detail in Sec. Ill is based on a 
modification of Stuart's model that takes into considera- 
tion both the orbital velocity and the steady-state 
cylindrical geometry. We briefly review the snowplow 
model and summarize the results for point releases. 

These results will be useful for comparing with the 
results of the venting model and will be used in Sec. V 
where we calculate the deposition of ionization from a 
knowledge of the time variation of the cloud size, R(t). 

A. Basic assumptions and equations 

The central assumptions of Stuart’s model are 
that a spherical cloud of gas, consisting of released 
gas plus swept-up ambient atmosphere, retains a self- 
similar or shape-preserving density profile, that the 
velocity of the gas is radial , and that the velocity 
increases linearly with radius from the center of the 
cloud. In the case of release at orbital velocity, the 
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above assumptions are modified. Ve assume that the 
expanding gas cloud remains spherical in shape and that 
its center moves with a velocity u(t). We also assume 
that the gas cloud completely picks up the ambient air 
by incorporating all of the atmosphere that crosses its 
boundary into the cloud. The cloud does not disturb the 
ambient atmosphere exterior to the cloud and no ambient 
molecules can pass through the cloud unaffected. The 
appropriate equations conserve the mass, momentum and 
energy of the expanding and moving gas cloud consisting 
of the released gas and swept -up ambient air. The 
qualitative results of the model are believed to be 
reliable during the early phase of the cloud's dynamics 
prior to times when diffusion is the dominant process. 

A detailed analysis resulting from these assumptions is 
described in Ref. 3. 

In the Stuart snowplow model and its extension 
described above, all lengths are found to scale with a 
characteristic length called the equal-mass-radius, a Q . 
This length is the radius of a spherical volume of ambient 
atmosphere with mass equal to the mass of the released 
gas, M q . In terms of the ambient atmospheric density, 

P a , the equal-mass-radius is defined by a Q • (3M 0 /4irp a )^. 
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T'lis model appears to be appropriate when 3 < a Q /X < 30 
where X is the ambient mean-free-path between collisions. 
At higher altitudes where a Q < 3X a kinetic treatment 
would be necessary. At lower altitudes where a Q > 30X, 
it is difficult to consider the released and ambient 
species as being well-mixed, and other gasdynamic 
phenomena neglected in the snowplow treatment , such as 
shock formation and wave generation, would be important. 

We summarize the equations that lead to a specifi- 
cation of the snowplow model. The similarity condition 
for the mass density, P(r, t), becomes 

p(r, t) - [M(t)/R 3 (t)] f(|) (1) 

where M(t) is the total mass of the cloud including 
swept-up ambient atmosphere and R(t) is the radius of 
the cloud. The dimensionless space variable t is given 
by 

t - [r - x(t)]/R(t) (2) 

where x(t) is the three-dimensional position of the 
center-of-mass of the cloud and r is any point in 
three-dimensional space. Equation (1) describes a 
gas cloud that retains its density profile while 
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expanding and translating. The shape function f(t) 
has the properties 


Jf(t)d 3 | - 1 ; 

Jtf(f)d 3 | - 0 . 


(3) 

(4) 


The second moment is a time- independent (but shape- 
dependent) constant a; 


a 




For a constant density profile, f(|) » 3/4* and a * 0.6. 
For a spherical shell with all of the mass concentrated 
at the radius R(t), f(£) * 6(|f| - l)/4ir and a * 1.0. 


The gas velocity at any point is assumed to be 
given by 


v(r, t) » u(t) + R(t)t(r, t) 


( 6 ) 


where 


u(t) 


dx 

dt 


(7) 


and the dot represents differentiation with respect to 
time. Thus the hydrodynamic derivative is given by 


D 

Et 



+ v 


7 - 



( 8 ) 
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because a point moving such that t is constant moves 
with the local "fluid” velocity. 

The conservation equation for mass is written 

§f ♦ P7 • v - D n V 2 p + S(r, t) (9) 

where S is a mass source funct ion due to sweeping up of 
the ambient atmospheric mass and D R is a diffusion 
coefficient. Snowplow models set D r ■ 0 in order not 
to violate the density similarity assumption. We will 
discuss the effect of diffusion on changing density 
profiles in Sec. II-C below. Due to the density similarity 
assumption, the mass source term must be given by 

S(r , t) - (M/M)p(r , t) (10) 

t where M is the rate at which ambient atmospheric mass 

is swept up as the cl ud expands and translates. For 
spherical clouds wit., ut translation 

t 

M - 4irR 2 p a R (11) 

With translation, a more complicated expression, 

' dependent upon model assumptions, is necessary. In 

general, with translation, the volume of swept-up 
atmosphere, V* - (M - M o )/p a> is greater than the 
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volume of the gas cloud, V ■ 4irR /3, 

The momentum equation Is written 

P(r, *) if + S(r, «) ▼ <*. t) - -7p(r, t) (12) 

where p is the pressure in the cloud and the mass 
source term accounts for the momentum of the picked-up 
atmosphere. The internal energy density of the gas 
is given by p/(y - 1) where y is the ratio of 
specific heats. The total internal energy of the cloud 
is 

U - f d 3 ? (13) 

A' 

where the integral is taken over the volume of the cloud. 
Following procedures similar to those applied by Stuart 
as discussed in Ref. 3, the equations governing the 
translation and radial expansion in the cloud can be 
brought to the forms 

g- (Mu) - 0 , (14) 

aR HE (MR) * 3(y - 1)U - 3p s V (15) 


20 



where a is the second moment of f and p. is the limit 
of the cloud pressure approaching the surface from within 
and has been assumed to be spherically symmetric. Note 
that Mu is the total momentum of the cloud and that both 
M and u are functions of time. 

Equation (15) shows that the expansion of the cloud 
is driven by (y - 1) times its internal energy, proportional 
to the cloud pressure, and is retarded by the pressure at 
the surface of the cloud. In our model, we take p e to 
be the atmospheric pressure p & . These equations represent 
a quantitative description of the basic snowplow model. 
Reference 3 contains discussions of the treatment of 
these equations, their relationship to other author's 
work, the effects of various treatments of the pressure 
at the surface of the cloud, extensions of the derivations 
outlined here to incorporate different values of the 
molecular mass and y for the released gas and atmospheric 
gas. 


We can replace the internal energy U, defined by 
Eq. (13), in Eq. (15) by using conservation of total 
energy , 

M - M p n 

E-U + K-E 0 + — p - -3 — ^ , (16) 

ft 
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where K is the kinetic energy given by 


ip 1 [ „£ 2.3 

K 2 / Pv d r 

J V 



( 17 ) 


E is the total initial energy of the cloud, (M - M )/p 
a V’ is the volume of the swept-up atmospheric gas, 
and p_/(y - 1) is the internal energy density of the 
ambient atmosphere. Because we are including the 
possibility of a nonzero translational Velocity for 
the cloud, the volume, V ' , of swept-up atmosphere may 
be greater than the volume, l/, of the cloud itself. 

The initial energy of the cloud, E , derives from 

o 

the translational velocity of the release vehicle and 
from the chemical and/or thermal energy of the exploding 
gas. The chemical energy, initially released as heat, 
is transformed into kinetic energy of expansion in a 
time scale which is short compared to that of the other 
processes in this model. Thus 


E 


o 



2 



• 2 

+ aR 0> 


( 18 ) 


where il is the initial released mass, u is the initial 
o o 

translational velocity, and R q is the initial expansion 
velocity. 
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Equation (15) becomes: 


aR ^ (MR) - 3 (y - 1) [e q - K + (V 1 - I/) • ( 19 > 

Ambient atmosphere has been swept up and compressed by 
the cloud and the term ( V' - V)p_/(Y - 1) is the portion 

ft 

of the swept -up atmospheric thermal energy that is 
available to drive the expansion of the cloud. If there 
is no translation, then V 1 = (M - M )/p * V, the ambient 

O CL 

atmosphere is not compressed and its internal energy 
does not contribute to the expansion of the cloud. At 
late times ( V - V )p„/(Y - 1) approaches a constant, 
thus contributing to Eq. (19) in the same manner as E Q . 

B. Results for point releases from rockets and Spacelab 

Equations (19) and (14) together with an equation 
for dl!/dt, such as Eq. (11), can be integrated numerically 
to produce the translational motion, expansion, and mass 
accretion of the gas cloud. Reference 3 contains the 
results of such integrations for various values of the 
input quantities E Q , y, a, and R Q . Here we show only 
typical results in order to illustrate the nature of the 
solutions and present analytic expressions that approximate 
the behavior of R(t) for R > a Q and have the correct asymptotic 
behavior as t ♦ 
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Figure 1 shows several examples of R(t) obtained 

by numerically integrating the equations. Both the 

radius and the time are scaled by ttye equal-mass-radius. 

The time scale has been drawn under the assumption that 
• 

R q * 1 km/s for all cases. The heavy solid curve is 
the nomin 1 case with y • 1.4 and a » 0.6 corresponding 
to a uniform density profile. The dashed curve shows 
the influence of the assumption of an extreme density 
profile, that of a spherical shell with all the mass 
concentrated at the outer edge characterized by o * 1. 

It does not expand as fast as the constant density 
profile because the available energy drives all particles 
at the same velocity. Likewise, the outer edge of a 
gas cloud with a density profile more peaked at the center 
with a < 0.6 would expand more rapidly than shown by the 
solid curve. 

The four thin solid curves show the effect of various 
values of y. As y decreases, the number of degrees of 
freedom increases and a larger fraction of the energy is 
stored as internal energy per particle and less is available 
for kinetic energy of expansion. The curve labeled y * • 
corresponds to the assumption that none of the available 
energy is stored as internal energy. Regardless of the 
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value of y associated with the released gas, as R >> 1 
the cloud radius approaches a curve corresponding to 
the value of y of the ambient atmosphere (as shown in 
Ref. 3) because the majority of the gas clouds' particles 
are atmospheric particles. 

The dotted curvj represents the original Stuart 
snowplow without the effect of translational velocity. 
This curve corresponds to a release from a rocket in 
which the translational velocity is small compared to 
the initial expansion velocity R q . All of the curves 
initially have the behavior R * R Q t . The snowplow 
without translation begins to decrease its rf*te of 
expansion as R a Q and the cloud has picked up an 
atmospheric mass comparable to the mass of released 
gas. 

The case of release at orbital velocity is qualita- 

tively different. Figure 1 shows that R increases from 

its initial value before R * a^ and reaches a maximum at 

around R - a Q before decreasing. The explanation for 

this behavior depends upon the fact that the velocity of 

the atmosphere relative to the shuttle corresponds to 

a large kinetic energy per particle amounting to 5.1 eV 

fox oxygen atoms. This energy input rapidly exceeds the 

*2 

original chemical energy, aM 0 R 0 / 2 > belore the swept-up 
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mass approaches IT or R ■* a ,, The result is an increase 

o o 

in the internal energy of the cloud leading to its rapid 
expansion. This rapid expansion leads to a rapid rate 
of sweeping up of atmospheric mass. The translational 
velocity of the gas cloud drops rapidly with a resultant 

i , 

significant reduction in the rate of energy input to the 
cloud. 

Figure 2 shows the typical behavior of a 100 kg 
gas release at orbital velocity at 200 km altitude. 

The time and length scales are given in seconds and 
kilometers, respectively, based on the derived value 
of 3.65 km for the equal-mass-radius. The dashed 
curve for the cloud radius is the same as the heavy 
solid curve in Fig. 1. The solid curve labeled tempera- 
ture corresponds to the average internal energy per 
particle with an assumed mass of 20 atomic masses. The 
curve shows the initially cold expanding gas rapidly 
heating due to the atmospheric energy input reaching a 
maximum of 5600°C by 2 s aft^r release. This rapid 
heating leads to the rapid expansion and subsequent 
cooling of the gas cloud. The temperature asymptotically 
approaches the ambient atmospheric temperature. The 
distance traveled Is,/.. * a\ * is shown by the 
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broken curve. The cloud at first coasts at orbital 
velocity but rapidly decelerates when R a Q reaching 
2/3 of its asymptotic travel distance within 3 s. Note 
that by 4.5 s after release the cloud diameter 
exceeds the distance traveled while the cloud is still 
expanding rapidly and its forward motion has nearly 
stopped. The fact that the cloud size rapidly exceeds 
the distance traveled suggests that the model assumption 
of a spherical cloud shape is justified because any 
asymmetry that may be introduced by the translational 
velocity rapidly becomes less important . 

Indeed, the fact that after a few seconds the 
gas release at orbital velocity appears to behave as 
a stationary point release of greater initial energy 
can be used to obtain an accurate approximation to 
R(t). In our model, when u(t) < R(t) or MR > M q u 0 , we 

« O • 

take M « 4trR p R and then V - V ■ AV ceases to grow 

li 

and M(t ) may be written 

M(t) - M + 4 it p R 3 + Al/p o (20) 

O O St Or 

where the last term is a constant generally smaller 
than the first and much smaller than the second. The 
total energy available for driving the late-time 
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expansion may be written 


Et - f 0 * iV y-rr (21) 

where E q is given by Eq. (18). The last term can be 
shown to be smaller than the first by writing it as 
(3al//4*a 3 ) C^ 0 C^/y(y - 1)] where C ft is the atmospheric 
sound speed. 

Thus, for UR greater than M q u 0 , Eq. (19) can be 
written 


aR ^ (MR) - 3( y - 1) (E t - K) 


( 22 ) 


where E T is a constant given by Eq. (21) and K is the 
kinetic energy given by Eq. (17). As t * «*, u(t) 0. 

M(t) - M o R 3 /a^. and K * ictMR 2 . If we define a non- 
dimensional time by 


s * R Q t / a Q (23) 

• 

where R Q is the initial expansion velocity, we find 
that for large s , the solution of Eq. (22) has the 
asymptotic form 

0 4 

R * a Q Cs * as s ■* » (24) 


where 


2S 



(25) 


C « 


f 


25( Y - li 

4 Y 


1/5 


/ 2E T \ 

V« M oV 


1/5 


The first term has the numerical values 1.123 and 1.201 
for y * 1.4 and 5/3, respectively. If we neglect the 
small term proportional to AV in Eq. (21), then the 
coefficient C can be written 



(26) 


By direct comparison with the numerical integration of 
Eq. (19) from s * 0 , we have found that an analytic 
expression of the form 


R(t ) » a Q C 




(27) 


provides an excellent approximation to R(t) for R » a 

o 

with C defined by Eq. (26) for both the stationary ana 
orbital release cases with u Q ■ 0 and 7.8 km/s, 
respectively. We shall make use of the expression in 
Eq. (27) for R(t) later in this report. 
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C. Snowplow model with diffusion 

At late-time, gas clouds released in the atmosphere 
have a Gaussian density profile and expand by diffusion. 

In this section we explore how rapidly the density profile 
of a spherical point release can change in profile from 
the self-similar profile assumed in the snowplow model 
to the Gaussian profile at late time. We find that this 
change can take place surprisingly rapidly. 

In this work we modify the snowplow model to include 
the fact that at late time the gas cloud can be modeled 
by a Gaussian density profile diffusing into the ambient 
atmosphere. The r „lt is a model which initially 
follows a snowplow expansion and evolves into a Gaussian 
diffusing cloud. Obviously, we must relax one of the 
central assumptions of the snowplow model, specifically 
that the cloud retains a shape-preserving density profile. 
We will treat the effect of diffusion as if it were 
superimposed on the expansion process modeled by the 
self-similar snowplow expansion. In a self-similar 
expansion the size of the cloud is given by R(t), a 
time-dependent scale length of the cloud. If the self- 
similar shape were a spherical constant-density profile, 
for example, we would take R(t> to be equal to the radius 
of the sphere. 
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We use this same R(t) to scale the continuity 
equation for mass, Eq. (9), but now we retain the diffusion 
term. Ve restrict this discussion to the stationary 
point release case for which f(r, t) ■ r/R(t), the fluid 
velocity of the gas is v ■ rR/R, and R(t) is the charac- 
teristic scale length of the expanding cloud if there 
were no diffusion. The last term on the right-hand side 
of Eq. (9) includes the effect of sweeping up ambient 
mass according to Stuart's original model and is given 
by Eq. (10). 

By making a suitable change of space and time variables 
we can define a function f[t(r, t), x D (t)] which is analogous 
to the shape function of the self-similar model. If we 
define t d by 


t D 



4D dt' 
n 

R 2 (f ) 


(28) 


implying that diffusion begins at t ■ t D , then the 
density can be written 


P(r, t) 


[s4(t)/R 3 (t)] 



t), T D ( t )] 


(29) 


With these substitutions Eq. (9) can be simply written 
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A 


( 29 ) 


ii_ » l v 2 f 

3t d 4 


where we have used the following: 


$ • v « 3ft/R ; 


3p «I f 

3t r 3 


,4 1 „3 St V „3 3t ' 


VP-^7 C 1 . 


JL v 2 f • 
R S E ' 


- - !tr/n ; 


3t D 4D n 

9t R 2 (t ) 


Equation (29) is the diffusion equation for f(t, t), 


which has the solution 


. t d ) - J a 3 r 


. exp t a - t,)2/T t>] 


t(Z\ 0) 


(ttt d )' 


(30) 


In terms of this function we may write the density at 
any time t > t n as 
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where Rp ■ R(tp) and Mp * M(t D ). For Tjj << 1 we have 
the limit of little diffusion and the density reduces to 
that given by the self-similar model, 

ofr - JL ^ o (r ^ t ) 

p(r * Mp p \ r R * t D/ 

or p(r, t) ■ f (r/R, 0) 

R J 

At long times after diffusion has been turned on, the 
similarity function f approaches a Gaussian distribution 
with an expanding radius and constant volume: 

f(5, t d ) -► (irr D )” 3 / 2 exp(-5 2 /T D ) (31) 


In order to illustrate the effect that diffusion in 
£ has on an assumed density distribution, we have applied 
Eq. (30) to two specific examples. The first is a uniform 
density out to a finite expanding radius; we refer to 
this profile as a flat-top distribution. The secoud is 
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an expanding spherical shell with all the mass concentrated 
at the surface of the shell. For the flat-top distribution 

UK, 0 ) = 3 / 4 tt ; £ < 1 

fU, 0) = 0 ; £ > 1 ; 

and for the spherical shell 

fU, 0) » <5(£ - l)/4ir . 

The normalization has been chosen so that each contains a 
unit volume. 

We have found that the profile approaches a Gaussian 

profile more rapidly if the Gaussian with which it is 

1/2 

compared has a Gaussian radius of the form (8 + x^) ' , i.e., 

UK, x D ) - [ tt (6 + x D )] -3/2 exp[- £ 2 / ( 8 + x D )] . (32) 

In order to illustrate the change of the flat-top and 
spherical-shell distributions toward the Gaussian form 
given in Eq. (32) at various values of x D , we have 
renormalized both the radius and density by defining 

z - £/ tO . 693 (8 + t d ) 1* , 
g(z, x D ) - [ir(8 + t d )] 3/2 f(£, x D ) . (33) 
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The limiting value of g(z, t d ) for large t d is 
2 

exp(- 0.693 z ). The scale factor 0.693 is chosen 
so that g(l, «) is 0.5 

Figures 3 and 4 show g(z, t^) for various values 
of Tq for the flat-top and spherical shell initial 
distributions respectively. The value of 6 in the 
scalings given by Eq. (33) are different for the two cases. 
It is appropriate that 8 be unequal to zero because in 
both cases the initial shapes have finite spatial extent. 

We have chosen 8 so that g(l, 1' is approximately equal 
to one half; for the flat-top 8 * 0.5 and for the spherical 
shell 8 ■ 0.75. The larger value of 8 for the spherical 
shell results from the fact that all of the mass is 
initially at £ ■ 1 for the shell while the mass is 
distributed over 5 1 1 for the flat-top. 

The vertical dashed lines in Figures 3 and 4 indicate 
the position of the initial radius in the newly-scaled 
radial coordinate z. For r D greater than one, the scaled 
position of the initial radius continues to move to the 
left and approaches z - 1.2 as + « indicating the 
increase in size of the cloud due to diffusion as 
shown by the ■ * curve. For the flat-top in 
Figure 3, the horizontal dashed lines indicate the 
height of the original flat-top in the newly-scaled 
density with amplitude [ir($ + t D |p^ 2 f(5, 0). 
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Examination of Figures 3 and 4 indicates that the 
density distributions are approaching Gaussian-like 
shapes for x^ > 0.2 and 0.75 for the two cases, respectively, 
and have essentially the asymptotic shape given in Eq. (33) 
for x^ > 0.4 and 1.0 respectively. We will discuss 
specific evaluations of x^ below. 

The approximation of f(C, x D ) by its limiting 
form [ir(0 + x D )]“* * exp[- £ 2 /(B + x D >] is equivalent to 
approximating the density distribution of the cloud by 
a Gaussian density profile 

p(r, t) * M ( t ) [na 2 (t)] _3/2 exp [- r 2 /a 2 (t)] (34) 

with a time-dependent Gaussifn radius, a(t), given by 

a 2 (t) - R 2 (t) [fi + r D (t )] (35) 

where x D is given by Eq. (28). Although the actual shape 
at early time is not Gaussian, in the case of the flat- 
top distribution we find that the distance at which the 
density is 1/e of its value at r * 0 is within a few 
percent of R(t) for x D (t) <0.5 and of course is 
indistinguishable from a(t) given by Eq. (35) for 


x D > 0.5. 



In the preceding, we have discussed a procedure for 
introducing diffusion within the framework of a snowplow 
model. The effect of diffusion is to transform an initial 
self-similar profile for the density into a Gaussian 
profile. The treatment so far does not depend on the 
the time-dependence of R(t) which we take to be determined 
by the snowplow model without diffusion. 

With a normalized time defined by Eq. (23), we 
define a function $(s) such that 

R(t) - a o 4>(s) . (36) 

With u » 0 we have AV » V* - V ■ 0 and the original 

Stuart snowplow equation has the form of Eq. (22) with 
• 2 *2 

E t ■ aM 0 R o /2 and K ■ oMR /2. Our normalizations for 
0 and s are different from those of Stuart; with 
our normalizations <J> satisfies the differential 
equation 

♦ as [a + $ 3 ) U] - [i - (i + <i> 3 ) (dd»/ds) 2 ]. 

The solution of this equation has the following 
properties for large and small s: 
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$(s) **“ s 


for s << 1 


<J> ( s ) - |e.25 1 j * s 0,4 for s » 1 . 

At early time the cloud expands with the initial expansion 

velocity R q while at later time the cloud radius is pre- 

CD 4 

dieted to expand as t * . At very late time when the 

snowplow description is no longer appropriate, the cloud 

expands as t®’® due to diffusion. It is apparent that 

once the appropriate value of y of the released gas is 

specified, the snowplow model for the expansion of a 

gas cloud determines the time-dependence of the radius 

R(t) of the cloud in terms of the two key parameters a Q 

and R . 
o 

Once the functional dependence of $(s) is known, 
then we may define a normalized diffusion function t(s) 
by 


x(s) 



(37) 


The diffusion parameter t d defined in Eq. (28) can be 
expressed in terms of t(s) as 
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( 38 ) 


4D n r_ _ 1 

T D (t) - l T ( s > - t ^ s D^J 

O O 

where s ■ R t/a and the diffusion is assumed to commence 

at t ■ t D . For a gas with y * 5/3 which is appropriate 

for barium atom vapor, we show the resulting normalized 

diffusion function 7(s) in Fig. 5. The solid curve is 

the correct value of x(s) only as long as <t>(s) is given 

by the snowplow model. If the snowplow model were to 

o 

be turned off at a specific value of s, then <j> in the 
denominator would be a constant and <J> would increase 
linearly with s. Examples are shown by the dashed curves 
in Fig. 5 if the snowplow were stopped at s * 18 and 
46. 

As an example of the time scale in which the profile 
can change from a flat-top distribution to a substantially 
Gaussian shape when ■ 0.2, we use parameters associated 
with a 48 kg barium release (called Spruce) at 190 km 

7 

altitude. For this cloud, detailed modeling has shown 

that D_ ■ 0.052 km^/s, a^ ■ 1.26 km, and R„ « 1.15 km/s 
n o o 

so that the numerical value of the coefficient of the 
brackets in Eq. (38) is 0.144. Thus t d has the value 0.2 
when t(s) - ^( s j)) " 1»4. Figure 5 shows that if we take 
s D as small as 1, could reach 0.2 as rapidly as 
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* 7 

t * sa o /^ 0 * 4.4 s with s - 4. Our detailed analysis 

of the observed profile of the neutral cloud has indicated 

that it is in substantial agreement with the scpIo size 

and profile obtained by adding diffusion to a snowplow 

description of the expansion in radius of the neutral 

clotid. 
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III. VENTING MODEL AND EQUATIONS 


In this section we describe a quantitative model 
that conserves the mass, momentum, and energy of the 
expanding and flowing gas in the wake downstream of a 
canister traveling at orbital velocity and venting 
gas. In the frame of reference moving with the canister, 
the flow is in steady-state. In accordance with snowplow 
models, we assume that all of the ambient air that 
encounters the wake and crosses its boundary becomes 
incorporated into it. The wake does not discurb the 
ambient atmosphere exterior to it. The quantitative 
results of this model are believed to be reliable beyond 
one wake radius downstream of the canister and prior 
to times when diffusion is the dominant process. The 
reasonableness and implications of these assumptions are 
discussed below. 

In the point release problem treated previously 
we found that all lengths scale with a characteristic 
length, the equal-mass radius, a Q . Likewise, in the 
steady-state wake problem all lengths scale with a 
different characteristic length a. 
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This length is defined as the radius of a cylindrical 
volume of ambient atmosphere that contains a mass per unit 
axial length equal to M Q /V a , i.e., 

a " <Vir y a P a >* * (39) 

* 

where M is the rate at which mass is being released 
o 

from the source. It is also the radius of a disk moving 

at the source (shuttle) velocity, V , in the ambient 

atmosphere through which the ambient mass flux is equal 

to M . 
o 

As in the previous section, we take this venting 
model to be most appropriate in an intermediate altitude 
regime defined by 3 < a/X < 30 where X is the ambient 
mean-free-path between collisions. At higher and lower 
altitudes the assumptions of the model are less justified 
and require different approaches. 

Figure 6 shows the altitude variation of the number 
density, n, mass density, p , and mean free path, X, of 

fit 

the ambient atmosphere based on the CIRA 1972 model 

atmosphere. The solid curves correspond to a moderate 

exospheric temperature T ■ 1200°. The mean-free-path 

\ ■ (na)” 1 is derived from an assumed typical collision 

- 19 2 

cross section c ■ 5 x 10 m independent of species 



or temperature. The variation of the atmospheric density 
for more extreme values of the exospheric temperature, 700° 
and 2000°, is shown by the dashed curves. 

Figure 7 shows contours of constant characteristic 
length a in the altitude versus mass- venting-rate plane. 

• X 

Because the horizontal axis is logarithmic and a * M Q a , 

all of the solid curves for T „ ■ 1200° have the same 

ex 

shape and are displaced horizontally from each other by 
an appropriate amount. The two dashed curves drawn for 
a ■ 1 km indicate the shapes of the contours for T ■ 700° 

v A 

and 1200°, respectively. These curves can be used to 
estimate the value of a for different model atmospheres. 

Figure 8 shows contours of constant values of a/X 
in the same altitude versus mass-venting-rate plane. 

Again, the solid curves all have the same shape and 
represent appropriate horizontal displacements. We 
note that the intermediate altitude regime, 3 < a/X < 30, 
covers a broad range of mass-rates and altitudes of 
interest for the experiments proposed for Spacelab. In 
this intermediate regime, ambient atmosphere is picked 
up and mixed with the gas cloud, and neither the high- 
altitude kinetic approach nor the low-altitude continuum 
fluid dynamics approach is valid. 
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Following the approach of our previous work 
outlined in the previous section, we have further 

5 

modified the treatment of Stuart and assumed a 
self-similar expansion of the cloud as a wake behind 
the source. We assume that the density of the wake 
is given by 

P(r. x*- f (£)M(x)/R 2 (x)V(x) (40) 

where x is the axial distance in the direction 
downstream of the source, r is the radial distance from 
the x-axis, R(x) is the radial extent of the wake at 
a distance x downstream, 5 * r/R(x), f(S) is the 
dimensionless self-similar radial density profile, M 
is the mass flux in the wake through the disc of 
radius R(x) centered on the x-axis, and V(x) is the 
x-component of the fluid velocity at axial distance x. 

As in the point release problem, an ad hoc assumption 
is required about the fluid velocity flow field. We 
have assumed that the x-component of the velocity, V(x) 
is independent of radial distance and that the radial 
velocity component is given by 

V r • v <*> KTJT a! ■ < 41 > 
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Equation (41) insures that the velocity at the surface 
of the steady-state wake is parallel to the surface. 

This model assumes that the only interaction 
between the gas cloud wake and the ambient atmosphere 
is the sweeping up of all the ambient atmosphere encountered 
by the cloud. Shocks and diversion of fluid flow in the 
surrounding atmosphere external to the cloud are neglected 
in this model. We also neglect the effects of ionization 
and condensation within the wake. 

The imposition of contraints on the velocity as 
in Eq. (41) and on the density profile in Eq. (40) is 
somewhat artificial as they have been imposed to make 
the problem mathematically tractable. These constraints 
lead to some inconsistencies between the model and the 
actual wake of a steady-state release. As with all 
mathematical models, it is hoped that the necessary 
inconsistencies are small enough that they do not 
destroy the value of the model. In this case, the 
model assumption about the velocity in Eq. (41) 
obviously cannot be satisfied near the point of 
release for nearly all modes of gas release. In fact, 
as we shall see later, Eq. (41) cannot be valid near 
the point of release for any configuration because of 
conservation of 'ass, momentum, and energy. 


45 


The equations for local conservation of mass and 


momentum are 


(P) + $ • (P$) « S(x, r) (42) 

and 

(P$) + $ • (P^) + h = t(x, r) . (43) 

In Eq. (42), S(x, r) is the mass source function due to 
the sweeping of atmospheric mass and the mass injection 
from the shuttle source. In Eq. (43), p is the gas 
pressure within the wake and t(x, r) is the momentum 
source due to the swept-up ambient atmosphere having 
velocity in the x-direction in the rest frame of 
the shuttle. 

The velocity vector defined above, is the 
fluid flow velocity in the rest frame of the shuttle 
and is not the velocity of the released gas with 
respect to the atmosphere as in the previous section. 

The corresponding conservation equation for energy is 

(u + | pV 2 ; + $ • [$(u + | pV 2 + p)l - Q(x, r) (44) 


where 
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u 


Y - 1 


(45) 


is the internal energy density, y is the ratio of specific 

heats, and Q(x, r) is the energy source term due to the 

swept-up ambient mass with velocity V and the energy 

81 

injected with the mass from the shuttle source. 

Because we are solving for a steady-state solution, 
the time derivatives do not contribute in Eqs. (42), (43), 
and (44). The global equations for conservation of mass, 
x-component of momentum, and energy are obtained by 
integrating these equations over the volume of the wake 
envelope bounded downstream by a plane of constant x 
as shown by the dotted line in the sketch below. The 
result is 
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1 


r-4- ■ 


M(x) “ M q + V a P a TrR 2 (x) = M Q (1 + R 2 /a 2 ) , (46) 


M(x)V(x) - p a V 2 irR 2 (x) - *R 2 p a + 


[* 

I P(x, 

J A 


r)2irrdr - 0 ,(47) 


-< u a + k v a + Pa>V R2 + 


r R 

V(x) J (u + §p |$| 2 + p) 


2irrdr 


(48) 


where u and p. are the internal energy density and 

ct a 

pressure of the ambient atmosphere and E q is the rate 

at which the source is injecting energy along with the 

mass injection M . 

o 

The ambient temperature T is related to p. and 

cl cl 

u a b ? 


kT 


a 


a 


m 


a 


(Y a - 1) U a 


(49) 


where k is Boltzman's constant and in is the mass 

a 

of the ambient particles. Thus, u n + p in Eq. (48) 

a a 

can be written 


u « + Po 

a a 


a 


kT 


a 


Y - 1 m 

'a a 


a 


r 
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I 


By using the definition of the ch&racteristic length, 
a ■ M 0 /irp a V a , Eq. (45) to relate p and u, Eqs. (40) 
and (41) for the radial dependence of p and and 
defining 


dU 

dx 



r) 2Trrdr, 


(50) 


Eqs. (47) and (48) can be written 


• R 2 / kT a \ 

MV - M ~ (V + I 

0 a 2 \ a Va / 


dU 


+ (Y - i) ^ * 0 , 


(51) 


" M ° ? U + Vi ST ) * YV § + 5 «< v2 + “R 2 ) - v < 52 > 


2 /V" 


Y a k T ft 


where 


a 



f(SK 2 2ir£d5 


(53) 


is a dimensionless invariant determined by the self- 
similar radial density profile and 


R - V « 
dx 


(54) 
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The radial component of the momentum conservation 

Eq. (43) is treated analogously to the virial treatment 
5 

of Stuart. By taking the dot product with r on both 
sides of Eq. (43), where r is the component of position 
perpendicular to the x-axis, integrating over the volume 
shown in the previous sketch, and differentiating with 
respect to x, one obtains. 

aR A (MR) * 2 (Y - 1) - irR 2 p a (55) 

The set of Eqs. (46), (51), (52), and (55) 
determines the four quantities R, M, V, and dU/dx and 
constitute the venting model. 
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IV. 


RESULTS OF THE VENTING MODEL 


In this section we illustrate the properties of 

the solution of the equations derived in the previous 

section. We restrict our discussion to the special case 

Y ■ y « Y and will not distinguish between the mass of 
a. 8 

the gas species, m , and the mass of ambient species, m , 

& a 

by setting them both equal to 20 atomic mass units as 

representative of a mean molecular weight in the thermosphere. 

The choice of the molecular mass has no impact on the 

dynamics when y = y but enters only when determining 
a g 

an average "temperature" defined as 

m V , TT 

T - (lf - w i ■ (56: 

A. Analytic characteristics of the solution 

In order to examine the nature of the solution of 
the venting model equations, we first note that the ambient 
energy density and internal energy appear in the same 
combination in Eqs. (51) and (55). Hence, we define 

U' - (Y-l) g - «R 2 P a . (57 

By using this definition of U' to eliminate dU/dx from 
Eq. (52) and rearranging terms, the energy equation may 



be written 


M q kT a /p\2 

Y-l V a nT (a) (V a“ V) 


+ 5 [ S ' o V a (|) 2 - - 5 aiti2 - (58 > 

As R+®, V+V and the second term on the RHS approaches a 

£L 

constant, contributing to the energy balance in the same 

manner as E Q . This terms arises from the compression of 

the ambient atmosphere as it is swept up and serves as a 

source of energy to drive the radial expansion. Each 

2 

of the terms in the square brackets increases as R as 
R-+-® aud the leading terms exactly cancel leaving a constant 
as the dominant contribution of this term to the energy 
equation in this limit. 

The structure of Eqs. (51) and (55) is considerably 
simplified if we define new normalized independent and 
dependent variables 


2 2 
y = R/a , 

(59) 

1 

\ 

q = MV/M 0 V a ; 

(60) 

1 

X = aM 2 A 2 /.\:£ ; 

(61) 

L 

* 

! 

j 

A * U’/M 0 V a • 

(62) 

) 
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Equations (51) and (55) can now be written 


q - y - A, (63) 

y t^* 2qA ’ (64) 


respectively. Unfortunately, the energy equation does not 
simplify greatly. However, by multiplying Eq. (53) by M, 
using Eq. (46), defining the constants c and b, 


c 



7 — 2 
M V 
o a 


b = 


kT 

a 

m V 2 
a a 


(65) 


and using Eq. (63) appropriately, the energy equation (58) 

2 2 

can be written uniquely with no terms of the form q , y , 
or qy; 


q* * c + (l+c)y + ^ y(l+A) + A 2 - X . (66) 

The three equations (63), (64), and (66) determine 
the three functions q, A, and x as functions of y (and thus R). 
The equations in this form do not contain R explicitly and 
the dynamics do not depend explicitly on the value of a. 

In terms of these normalized functions, the axial velocity, 

V, and temperature, T, defined by Eq. (56), are given by 
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(67) 


V = V = 

a 1+y 


V 

a y+1 


?( 

M 


m. T J' 


iJo/RV 5 T \ = !!a v 2 (y-A)(by+A) , 

V- a / T a/ - k v a (y+1) 2 ’ 


( 68 ) 


where the last form is an explicit function of y once a 
solution A(y) of the equations has been obtained. 

One consequence of the venting model and the above 
equations is that solutions do not exist below a minimum 
value of R which depends on the various parameters. 
Accordingly, solutions can only be obtained for finite 
initial values of y. For any given set of gas release 
parameters and a sufficiently large value of y (or R), the 
model allows a family of solutions corresponding to different 
values of P. (or xi see Eq* (61)). The range of allowed 
values of R (or x) for a given value of R (or y) is 
determined by the energy equation. The maximum permitted 
value of x is given by Eq. (66) with A=0, or 


X = c + (1+ c + — j) y . 

On the other hand, if Eq. (63) is used to eliminate q 
fron Eq. (66), we obtain another form of the energy 
equat ion 
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(69) 


I±1 *2 . gX^bi yA + c +( i +c + |l_. b) y . x = o 

which is quadratic in L. The requirement that A be real 
puts a lower bound on the permitted value of v in order that 
the discriminant be positive ior sufficiently small y. 

The model equations have been integrated numerically 
for a variety of different allowed initial conditions. All 
of these families of solutions asymptotically approach the 
same behavior at large radius far downstream in the wake 
where the model assumptions are more valid. Some of the 
allowed values of initial conditions at small R are within 
reasonable ranges and some are not, but all such solution 
curves are less valid at small R because the model assumptions 
are less justified. 

Before illustrating the numerical solutions we 
obtain the asymptotic form of the solutions and give 
an analytic expression for the expansion velocity that is 
an excellent approximation to a particular numerical 
solution. By direct substitution into Eqs. (63), (64), and 
(66) or (69), it can be verified that the following forms 
represent asymptotic solutions of tue equations for the first 
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two terms in an expansion valid as y-*"» 


‘•‘o*?’ °(^) ' < 70 ' 

x - 24 o y + 2i 0 B + o(i) ; (71) 

V ' * V a i 1 - TT^T? ) + °{p) ’ < 72) 

T * T a * S V a IT + 0 (■?) '■ (73) 

where < * £7^ (1 + c) 

A q = (< + b) ' (1 - b) . v 74 ) 

B * c/2A q + (^- b + l) A q /2 . (75) 


All solutions of the model equations have the 
asymptotic form given by Eqs. (70)- (73) for sufficiently 
large y. However, the analytic expressions obtained by 
keeping only the first two terms in the expression (71) 
and setting A ■ A Q in Eqs. (67) and (68) are indistinguishable 
througn three significant digits from a particular set of 
numerical solucions for y 2 1. By remembering the definition 
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for x» Eq. (61) , we can write an explicit expression for 

R, 


_ K 2 v (R 2 /a 2 + B) 1/2 
2 a R 2 /a 2 + 1 ' 


where 

K 2 /2 = (2A 0 /a) 1/2 . 


( 76 ) 


(77) 


The expression (76) for R(R) has an asymptotic expansion 
for R>>a, 


V al ^|] ♦ 0 (si)_ (78) 

identical through the first two terms to the asymptotic 
expansion of all the solutions. 

The parameter b, defined by Eq. (65), is a normalized 
ambient temperature and has the numerical value b ■ 0.00679 
for T ■ 1000°K and V * 7.8 km/s. Hence it makes only a 

£L cL 

very small change to the value of and B defined by Eqs. 

(74) and (75). However, the form of t Q given in Eq. (74) 
with the term (1-b) in the denominator indicates that this 
venting model cannot be used to describe the gas cloud 
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produced by slowly moving canisters. Indeed, it is probably 
inaccurate for velocities approaching th*- ambient thermal 
velocity. Thus this model must not be used to describe 
venting releases from sounding rockets which have velocities 
in this range. This breakdown in this model is associated 
with the compression of the atmosphere that results from 
our description of the atmospheric mass pick up. 

The parameter c, defined by Eq. (65), is related to 

the stored chemical or internal energy per unit mass. For 

reasonable values of E /£L. c << 1 and does not have a 

o' o 

major affect on the solutions. The detailed analyses of 

7 

barium releases from rockets indicated that the best value 

of the initial expansion velocity of the barium vapor with 

a flattop distribution is * 1.15 km/s. This value 

• o 5 

corresponds to E Q /M 0 * 0.6 R Q /2 « 4 x 10 J/kg. We use 

this same value for t 0 /& 0 in the case of a vented release 

and obtain c ■ 0.0131. With these values for b and c, the 

o 

quantities b Q , B, and K V ft /2 have the values tabulated in 
Table I for the indicated values of y. 

B. Evaluation of the solutions 

The equations for the venting model have been 
integrated numerically. The next twelve figures show 
examples of representative solutions for various gas 
release parameters and assumed initial conditions. In 
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TABLE I. Numerical values 


Para net er 

r - 1.1 

Y * 1.4 

Y - 5/3 

A o 

0.0532 

0.1523 

0.2108 

B 

0.1528 

0.123 

0.140 

\ V a (km/s) 

3.60 

6.09 

7.16 


order to illustrate the sensitivity of the solution to 
a change in the value of each of the various parameters, 
we have compared the solution curve produced by the 
changed value of a single parameter with a "nominal" 
case corresponding to a standard set of parameters. 

Table II lists the values of the parameters for which 
representative solutions are shown. The first column 
of values represent the values associated with the 
nominal solution. Solutions associated with the other 
values are illustrated in the figures. 

The parameter 0 o is related to the stored 
chemical or internal energy and is given by C Q ■ /c V ft . 
It is the velocity at which all the mass would travel 
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TABLE II. Venting model parameter values. 


Parameter 

Nominal 

Value 

Other 

Values 

C Q (km/s) 

0.89 

3.09 

a 

0.5 

1.0 

FL (km/s) 

3.23 

2.23, 4.23, 1.93, 3.76 

Y 

1.4 

1.1, 5/3 

T fi (°K) 

1000 

0 


if all of the stored initial internal energy were 

converted into kinetic energy. The nominal value 

is appropriate for the specific energy associated with 

5 

barium vapor of 4 x 10 J/kg. The larger value of C Q 

corresponds to a value of specific energy twelve times 

larger than for barium. If applied to hydrogen gas, 

it corresponds to the energy content of gaseous H at 

* 

770°K and is 8.5% greater than the mean thermal velocity. 

The parameter a is the density distribution 
parameter defined for this cylindrical geometry case 
by Eq. (53). The values 0.5 and 1.0 correspond to 
a flattop distribution and cylindrical-shell 
distribution, respectively. (In the point release 
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case, a has the value 0.6 for a flattop distribution 
as discussed in Section II.). 

In Fig. 9 we show the nominal solution for ft(R) 
by a heavy curve and the solutions corresponding to 
the alternate values of C Q and a by lighter solid lines. 
In this and the subsequent eleven figures, the heavy 
curve always corresponds to the nomiual case. The 
other curves represent solutions corresponding to 
changes in only one or two of the parameters. These 
curves are labeled by the non-nominal value of the 
parameter. At the top of the figure the nominal values 
of these parameters are given . The values of all 
parameters not shown explicitly are those given as 
the nominal values in Table II. Each figure thus 
emphasizes the sensitivity of a solution to the change 
in a particular parameter. 

Returning to Fig. 9, the nominal curve is 
specified by the choice R Q • 3.23 km/s at R ■ R ■ a. 

In this and the next eleven figures the scale length 
a, defined by Eq. (39), is represented by A. This 
value for is the value produced by the expression 
(73) evaluated at R ■ a. The analytic expression (76) 
is indistinguishable from the numerically-integrated 


61 



solid curve. The two other curves behave like expression 

(78) for R >> a (in practice. R £ 2a) with the appropriate 

evaluations for K and B given by Eqs. (77) and (75). For 

C Q * 3.09 km/s the expansion velocity is only slightly 

greater for a given value of R/a due to the small additional 

amount of initial energy available for expansion. Note 

that there is maximum in ft(R) for R > a. If a larger 

initial value for ft at R_ = a had been chosen there 

o o 

would have been no maximum for R > a. We will illustrate the 
effect of the choice of ft on the solution in the next figure. 
For a mass distribution more peaked towards the outer 
edge, the expansion velocity approaches that given by the 
a = 1.0 curve corresponding to a cylindrical shell. 

The expansion velocity is less in this case because the 
energy available for expansion is distributed to all of 
the particles. The expansion velocity of the outer edge 
corresponding to a density distribution more peaked in 
the center with o < 0.5 would be greater than that given 
by the nominal curve, varying as a~^. 

Figure 10 shows solutions resulting from different 
choices for the initial value of ft Q . The value? illustrated 
correspond to the nominal value tl km/s. Both curves 
quickly approach the nominal solution curve which is well 
approximated by Eq. (76). For R > n :\, the solution is 
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fairly insensitive to the particular choice of R Q for 
reasonable values of this parameter. The asymptotic 
expansions given in Eqs. (70)-(73) and (78) depend only 
on the released gas parameters and are independent of 
ft . For the ft A ■ 2.23 km/s case, a smaller fraction 
of the total energy is contained in the kinetic energy 
term and a larger fraction is contained in internal 
energy with the result that initially the gas has a 
higher temperature, as we shall see below. This 
excess thermal energy rapidly expands the radius until 
the expansion velocity approaches the nominal solution. 

Figure 11 illustrates the dependence of R(R) 

on the assumed ratio of specific heats, > , according 

to the venting model equations. The two dashed curves 

represent solutions appropriate for the release of 

monatomic (y ■ 5/3) and polyatomic (y * 1.1) gases. 

In general , as y decreases, the number of degrees of 

freedom increases with the result that more energy 

is stored as internal energy and less is available 

for kinetic energy of expansion. The initial value 

of ft ■ 1.93 km/s for the y ■ 1.1 case is the value 
o 

obtained from Eq. (76) evaluated at R ■ a. The thin 
solid line with R 0 * 1.93 km/s illustrates the affect 
of a change in y on the resulting expansion velocity. 
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Equation (76) has the value 3.82 km/s at R 3 a with 
Y * 5/3. This value is below the minimum permitted 
value or' R Q at R * a according to Eq. (69). Rather 
than choosing a numerical solution with a larger 
initial value for R Q , we have shown a numerical 
solution beginning at R Q = 1.1a with a permitted 
initial value for R q of 3.76 km/s given by Eq. (76) 
evaluated at R * 1.1a. With this choice , the two 
dashed curves are indistinguishable from the analytic 
expression (76). Regardless of the value of y associated 

W 

with the released gas, when R >> a the effective value 

of y of the gas mixture will tend towards y associated 

a 

with the ambient atmosphere as shown in Ref. 3 for the 
point release case. The appropriate value ranges between 
1.4 and 1.67 depending upon the altitude of release. 

The dependence of the temperature, T, defined by 
Eq. (56) in terms of the average internal energy per particle, 
on the gas release parameters is given by Eqs. (68) and (73). 
An accurate representation for the nominal case for R > a 
is given by Eq. (68) with A - A . Figure 12 shows T(R) for 
the nominal case and the dependence on different initial 
values of R q . A smaller initial kinetic energy of expansion 
results in a higher initial temperature. The temperature 
resulting from an assumed expansion into a cold atmosphere 
is shown by the curve labeled T * 0°K. The presence of 
a warm atmosphere has only a very small effect on the dynamics 

04 


r 



of the expansion due to the small value of the parameter b. 
The primary effect is to raise the temperature of the gas 
cloud mixed with ambient atmosphere by T for R > 2a as 
shown explicitly in Eq. (73). Figure 13 shows the effect 
of y on the temperature. The greater the number of degrees 
of freedom, the lower the temperature. 

Figure 14 shows the axial velocity of the gas plus 
air mixture in the frame of the canister. As the gas cloud 
encounters the ambient air, it comes to rest with respect 
to the atmosphere which streams by the canister at speed 
V a = 7.8 km/s. The analytic.. 1 dependence for the nominal 
case is given adequately by Eq. (67) with A * A q . Figure 14 
shows that the velocity is rather insensitive to the choice 
of R q . Equation (67) with A = A q indicates that the velocity 
of the wake with respect to the atmosphere varies as 


wake 


V - V~ 
a 


1 1 *0 

i + R^/a^ 


a 


for 


R > a 


(79) 


which is rather insensitive to gas release parameters because 

A is small, 
o 

The shape of the plume behind the canister is obtained 
by recognizing that dR/dx ■ £/V and integrating 
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V(R’) 

R(R’) 


dR' 


(80) 


x 


/ 


By using the approximate expressions (67) with A = l ar.d 
(76), we obtain the explicit relationship 


x - x 0 “ 3 ( R ( R2 / a2 + B ) 


- (2A Q + B) a 2-n 


[R/a + (R 2 /a 2 + B)^j) 


(81) 


where x q is an arbitrary constant of integration. Asymptotically, 
for R >> a, the radius of the plume is found to vary as 


R(x) 


^(axV 


1 + (2A + 


Ln_(R/a_2 
DJ 2 9 

R /a^ 



(82) 


with K given by Eq. (77). Far in the wake x~ V t and 

a, 

0 5 

the radius of the plume increases as t ' . This asymptotic time 

0 4 

dependence is different from the t ' dependence f or the 
spherical point release case because of the different geometiv. 

Figure 15 shows the normalized plume radius as a 
function of normalized distance downstream of the canister 
for various values of the gas release parameters. The 
constant of integration has been chosen in each case so that 
R = a (* 1.1a for y * 1-1) at x ® a. The half-power 
dependence of radius on distance is obvious for R > 2a; 
the value of the y-axis intercept of the extrapolation of 
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the straight-line portion of the curves is equal to K. 
Figure 16 shows the same set of curves on a linear scale. 
Note that the vertical axis is expanded 9.3 times the 
horizontal axis. Figure 17 illustrates the sensitivity 
(or insensitivity) of the plume radius to assumed initial 
conditions . 

We expect that the wake model provides an adequate 
description of the characteristics of the vented gas once 
R > a. We expect that this criterion occurs wi .in a 
distance of order a downstream of the canister. We 
have no model for the transition region from the vicinity 
of the canister to the cylindrical region described by the 
wake model. Figure 14 shows that the released gas has lost 
half of its orbital velocity by the time this wake model is 
applicable and is decelerating rapidly. The distance. L, 
that the released gas travels with respect to the atmosphere 
as it grows in radius from R q to R is 
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where we have made use of Eqs. (76) and (79) for the 

approximate expression. Because for large R we have 
0 5 0 5 

R ~ x ‘ and R ~ t ' , we find that this model predicts 
that the distance traveled is not finite but increases as 
£n t. This result is different from the spherical point 
release model result which has the gas cloud traveling a 
finite distance. 

Figures 18, 19, and 20 show the normalized distance 
traveled by a parcel of gas located a distance x downstream 
in the wake for various gas release parameters and initial 
conditions. In these figures. L is arbitrarily assigned 
the value a at x = a. The three cur s in Fig. 19, all 
drawn for the same gas release parameters, indicate that a 
slab of gas initially having R = a at x = a , travels 0.2a 
fu^tuer (less) if the initial expansion velocity is 1 km/s 
less (greater) than the nominal value. 

In summary, the wake model derived here produces 

quantitative estimates of the primary characteristics from 

a continuously venting gas release traveling at orbital 

velocity. All lengths scale with the characteristic scale 

length a. Far downstream in the wake the cylindrical plume 

0 o 

radius expands as t ‘ . the shape of the plume as a function 

0 5 

of distance x downstream is R ~ K(ax) ‘ , the axial velocity 
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2 2 

falls as a /R and the distance traveled increases indefinitely 
as In t . a position x ~ 100a downstream, the typical 
plume radius is ~ 12a and the distance traveled with respect 
to the ambient atmosphere since R = a is AL ~ 4a. Close to 
the canister when R ~ a, the expansion velocity falls rapidly 
from ~ 3 km/s and the temperatur ' rapidly cools from several 
thousand degrees to the ambient atmospheric temperature. 

These results are dependent upon the snowplow model 
used to derive the appropriate equations. The most suspect 
assumptions are those relating to the self-similar radial 
density distribution and the uniformity of the axial velocity 
v.’th radius. It would seem reasonable to expect that both 
the mass pick-up and momentum pick-up in the wake are probably 
more concentrated toward the outer edges than assumed by this 
model. Near the canister, the outer edges of the plume 
wuuld oe swept back more rapidly than the central portion, 
particularly in the lower altitude end of the intermediate 
altitude regime. This effect is partially modeled by the 
c* * 1.0 curve corresponding to the assumption that all of the 
released and picked-up mass is concentrated at the outer 
°dge. The result is a small reduction in the radius of 
the plume. Far downstream where the axial velocity with 
respect to the atmosphere is small, V wa ^ e << R, R >> a, 
and p — p , the results of this model for a flattop 
distribution, a s 0.5, should be more valid. Of course, 
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at very late time, the released gas in the plume expands 
by diffusion. Th'' radial distribution can be assumed to 
evolve from a flattop to a Gaussian shape during the snowplow 
expansion as outlined in Section II-C. 


70 



V. CONDUCTIVITY ENHANCEMENTS PRODUCED BY BARIUM RELEASES 

Barium releases from rockets have been conducted in 
the ionosphere worldwide primarily to measure ambient 
electric fields using small releases (< 4 kg) at high 
altitude (> 200 km) and to produce enhancements in the 
ionospheric electron concentration using large releases 
(16 kg to 320 kg) at lower altitudes (< 200 km). The 
properties of the larger releases have been studied 
extensively using a variety of radar, optical, propagation 
path, and rocket probe techniques. As a result, many 
characteristics of these releases, such as the time development 
of the scale size of the neutral cloud, peak electron 
concentration, and distribution of barium ions perpendicular 
and parallel to the magnetic field are well documented 
and the behavior is understood. 

These large barium releases have produced major 
enhancements in the ionospheric field-line-integrated 
Pedersen (FLIP) conductivity. The ability to create such 
large FLIP conductivities is of interest for a variety 
of proposed Spacelab exper ments. In this section we first 
briefly describe the process of creating ion clouds from 
released barium atom vapor. We then apply the results 
of Sections II and IV for the expansion of gas clouds 



released in both the explosive and venting modes to the 
calculation of ion deposition from releases at orbital 
velocity. The results of these calculations are estimates 
of the scale size and FLIP conductivity enhancements that 
can be produced by barium releases from rockets and 
Spacelab as a function of altitude of release and amount 
of released vapor. 

A. Ion deposition for point barium releases 

In Section II we summarized the results oi a snow- 
plow model that describes the expansion of a spherical gas 
cloud in an intermediate altitude regime released either 
at orbital or negligible velocity. A barium thermite 
release typically produces an expanding cloud of barium 
atom vapor that exits the canister in 0.1 - 0.2 seconds 
and thus this event can be called a Ant release. Our 
understanding of the process by which jarium ion clouds 
are produced from the expanding vapor cloud and the resulting 

p operties of the ion cloud comes in part from a detailed 
7 

analysis of the Spruce ion cloud resulting from a 48 kg 
barium release at 190 km altitude, and in part from observations 
of many barium releases ranging in size from 16 kg to 320 kg 
and released at various altitudes between 150 km and 255 km. 

The description of the ion deposition as well as the 
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quantitative values for important parameters for the process 
are taken from Ref. 7. Here we summarize this process and 
describe the modeling adopted to obtain quantitative 
estimates for other altitudes and chemical payload weights. 

As the barium atom vapor expands, it encounters 
the ambient atmosphere. The outer radius expands initially 
with velocity R Q . We assume that diffusion changes an 
initial flattop distribution into a Gaussian profile 
whose Gaussian scale length increases according to the 
snowplow model. The number of barium atoms is depleted 
with an exponential time constant, x., by two competing 
processes: oxidation by ambient oxygen molecules and 

photoiou ization from excited metastable states. The 
cross section for the absorption of solar radiation that 
excites ground-state barium atoms into the metastable 
states is large. Photoionization commences only after a 
time dealy, t^ when the barium cloud has expanded 
sufficiently so that the metastable exciting radiation 
is able to permeate the cloud. 

As the barium ions are created, they become tied 
to the magnetic field line at their position of creation 
and cease their motion transverse to the magnetic f <_ ild, 

5 . After the passage of several x^ time periods, the 
bulk of the ionization has been created. The number of 
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ions on a given magnetic field lino remains constant beyond 
this time and the distribution of this field-line content, 

N (r), with radius transverse to § is a bell-shaped curve. 
We represent the characteristic size of the ion cloud by 
a scale length, a A , which is the distance transverse to 

-V _1 

D in which N c has fallen to e of its central peak value. 

The process described above has been modeled in 
detail in Ref. 7 for the Spruce ion cloud. The results 
of that model are the determination of appropriate values 
of a number of the required parameters so that the model 
produced results consistent with the various measurements. 
Table III presents a tabulation of the values of these 
parameters . 

A popular barium thermite reaction that produces 
barium vapor is 

CBa + CuO ► BaO + Cu + (C - 1) Ba 

where ; Is a number typically between 1.7 and 2.5. For 
Spruce with C * 2.5, the 48 kg chemical payload weight 
produced 3.56 kg of barium atom vapor which corresponds 
to a vaporization efficiency of 7.4%. Quantitative 
determinations of this efficiency are dif fuclt to obtain 
to much better than a factor of 2 and typical reported 
values range from 5% to 12%. The model value reported 
here is consistent with measurements of the properties 
of the neutral and ion clouds. 


(84) 
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TABLE III. Parameters for the Spruce ion cloud. 



Parameter 

Value 

Units 

Inputs to Model 



1. 

Mass of barium vapor, M q , 

3.56 

(kg) 

2. 

Vaporization efficiency 

7.4% 


3. 

Fraction of atoms ionized, f, 

0.42 


4. 

Ionization time constant, , 

10.0 

(s) 

5. 

Total ion inventory 

6.5 

(* 10 24 ) 

6. 

Atmospheric density, p . 

SL 

0.45 

(kg/km 3 ) 

7. 

Neutral diffusion coefficient. D , 

n 

0.052 

(km^/s ) 

8. 

Equal-mass-radius, a Q , 

1.26 

(km) 

9. 

Initial expansion velocity. . 

1.15 

(km/s) 

10. 

Onset of photoionization, t^. 

1.9 

(s) 

Derived Quantities 



11. 

Neutral cloud radius, R, at t = t^ 

1 . 59 

(km) 

12. 

Neutral cloud radius, a, at t * t^ 




+ T. 
l 

3.8 

(km) 

13. 

Transverse scale size, a x , 

2.7 

(km) 

14. 

Field-line content. N , 

V 

2.4 

ao 17 m- 2 ) 

15. 

Peak ion density at 100 s 

1.7 

(10 13 nT 3 ) 

16. 

Ion cloud length at 100 s 

7.8 

(km) 

17. 

Initial Pedersen conductivity 

27 

(mho ) 





The fraction, f, of the atoms that are ionized is 
equal to the ratio of the photoionization rate, v, to the 
total loss rate of barium atoms. The photoionization 
time, 1/v, has been consistently reported by Haerendel 

g 

et al . of the Max-Planck Institute to be of order 19 s. 


Rosenberg et al. “ of AFCRL report a time closer to 30 s. 

We choose 1/v = 24 s as being within 20% of the correct 

value. The oxidation rate is proportional to the molecular 

oxygen concentration, kn n , and was estimated by Rosenberg 

u o 

9 

et al . from observations of a series of barium releases 

at different altitudes. Because n n is considerabj ' lower 

u 2 

in the 1972 CIRA atmosphere model than in the U.S. Standard 

Atmosphere used by Rosenberg et al. , we have had to use 

a relatively large value for the rate constant, k = 1.25 
**16 3 

* 10* m /s , in order to match their observed loss rates. 
The fraction ionized is 


f 


v 


v + kn 


1 

1 + 3 x io 
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and has the value for Spruce given as entry 3 in Table III. 
The time constant for the loss rate of barium atoms is 


t . 
l 


v + kn. 


f 

v 


24 


1 + 3 x l0" 15 n. 


( 86 ) 
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and has the value given as entry 4 for Spruce. The total 
number of oarium atoms is * M 0 ^ m Ba w ^ ere m Ba is the 
atomic weight of barium. The total ion inventory given 
as entry 5 for Spruce is fN^,. 

If we neglect transport of barium ions perpendicular 
to the magnetic field, the conservation equation for 
the ions can be written 



vn b (r,z,t) 


(87) 


where v. is the velocity in the z-direction assumed to 
12 

be along $, n b (r,z,t) is the concentration of barium 
atoms, and r is the radial coordinate perpendicular 
to Based on the discussion in Sec. II, we take 
the distribution of barium atoms to be an expanding 
Gaussian of the form 


2 2 2 

n b (r,z.t) * - — 2 J 2 exp [-( r^ + z*)/R M (t)] e 1 

[rrR 2 ( t ) ) 


( 88 ) 


where t* is the length of time that the barium atoms have 

been exposed to ionizing radiation. By integrating Eq. (87) 

along the magnetic field we obtain an equation fox' the time 

dependence of the radial distribution of the field-line 

content N' (r.t) * / n.dz 
c v 1 


3X c (r.t ) 

3t 


v v _ r * / - 

T r o o 1 i 

exp t-r“/R“ ( t ) ] e 


tR 2 ( t ) 


i t 


(89) 



If we start the ionization process at time, after release 

we can integrate Eq. (89) in order to determine the radial 
profile of the field-line content of the ion cloud, 


N c (r,t) = vN t 


^ exp [-r 2 /R 2 (t ' ) ] (t ’ 

J TTR 2 (t') 

'i 


t. )/T. 
i" 1 


dt 


( 90 ) 


Entries 6-9 in Table III are related to the neutral 
cloud expansion. The first two are atmospheric numbers for 
190 km altitude and determine the equal -mass -radius. As 
discussed previously, it is appropriate to delay the start 
of photoionization until the neutral cloud has expanded 
sufficiently so that the metastable exciting radiation can 
permeate the neutral cloud. For the Spruce ion cloud, this 
delay time was considered a free parameter and adjusted until 
the time-dependent characteristics of the ion cloud matched 
the available data as described in detail in Ref. 7. The 
radius of the Spruce neutral cloud at the delay time t^ was 
1.26 a Q ••’1.59 km as given in entry 11. 

If we define a as the size of the Gaussian radius 

of the neutral cloud at the time t = t. + t. after release. 

i l 

i.^., at one ionization time constant after ionization has 
been assumed to begin, we may write the field-line content 
as 


7S 



( 91 ) 


fN T 

N (r ,t ) * I(r,t) 
ira 


where 


I(r 


-■* /' 

z i 


exp [~r 2 /R 2 (t ’ )] 
R 2 ( t ' ) 


“ ( t ’ - t.)/T. 
e dt 


(92) 


I(r,t) represents a normalized field-line content. The 
normalization is chosen so that 




I (r , t ) 27rrdr = 


exp [-(t 



(93) 


which approaches one as t -*■ ® . 

The result of the integration in Eq. (92) combined 
with Eq. (91) is the field-line content as a function of 
radius transverse to the magnetic field. From the profile 
we determine that for Spruce the radius at which the content 
has a value e~* of its central peak value is a^ « 2.66 km. The 
value of 1(0,®) is 1.65 which when combined with Eq. (91) 
determines the peak field-line content through the center 
of the cloud as given by entry 14 in Table III. 

Fcr a point release at orbital velocit •, vk 
neglect the translation of the neutral cloud and >-se the 
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approximation that the ions are deposited as though the 
neutral cloud were stationary. As pointed out in Sec. II, 
the distance that the neutral cloud travels is ,‘inite and 
quickly becomes small compared to the rapidly expanding 
spherical cloud. Hence the same procedure outlined above 
can be applied for this case with an appropriate specifi- 
cation for R(t) as is discussed in Sec. V-C. 

B. Ion deposition for venting-mode barium thermite 

releases 

Thermite reactions have been successfully used 
in order to produce lithium atom vapor trails released 
in a venting mode from both rockets and satellites . ^ To 
date, there have been no releases of barium atom vapor 
via a thermite reaction in the venting mode. Thus, we 
have no data regarding the nature of the ionization 
enhancement that can be produced by such a release. In 
this section we calculate porperties of such an enhancement 
by adapting the ion-deposition procedures applicable to 
the Spruce ion cloud to the trail-like neutral-cloud expansion 
model resulting from a venting-mode release. 

In principle, the procedures are similar but 
the details of the calculations are different due to the 
cylindrical geometry and the motion of the neutral cloud 
in the axial, x' , direction. We assume that there is a 
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critical radius, R^ , to which the cloud must expand before 
photoionization and loss of barium atoms commences. (We 
will discuss the criteria for determining this critical 
value as a function of altitude anu mass-venting-rate, 

M q , in the next section.) In the rest-frame of the shuttle, 
the venting-mode model presents a steady-state shape for 
the cloud-wake, R(x') with the origin chosen so that R(0) 

= R i , in which the gas has an axial flow velocity V[R(x')] 

■ V’(x') which approaches V as x’ -*• «. In the (x,y,z) 

a 

rest-frames of toe atmosphere, the shuttle travels in the 
negative x-direction, which we shall assume to be perpen- 
dicular to S, with velocity V = -V . At time t = t . , we 

So. X 

choose the origin of the x-axis to be located at the 
position where R = R^ . Due to the coordinate transformation 

x' = x + V (t - t ) , 

a _l 

the cloud radius, R , in the rest frame of th? atmosphere 
is 


H c (x,t) » R(x ' ) i R [x + V a (t - t.)] . 

After the passage of many t i time periods, the 
bulk of the ionizat on has been deposited and is uniform 
in x. Henceforth, we will concentrate only on the x = 0 
position and examine the t'me dependence of the ion 
deposition on field lines as the expanding neutral trail 
sweeps by that position. The cloud radius at x * 0 becomes 


(94) 
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a function of time only; 


R c (t) = RIV a (t ' t i )] • (95) 

Note that in this steady -state problem, the time t^ 

represei ts a critical time such that at x = 0, R * R. 

0 1 

at t = t.; it does not represent a time delay since 
release as in the point release problem. 

During the ion deposition we model the barium 
cloud trail as having a Gaussian density profile in the 
radial direction with a Gaussian scale length R,(t) = R(x' ) 
given by our self-similar snowplow expansion for the 
reasons discus ?d in Sec. II-C. The barium atom concentration 
that enters Eq. (87) at x = 0 for t > t^ is now written 

Nrp 0 2 o 

n h (y,z,t) = p exp [-( y “ + z )/R^]e (96) 

D VttR^ c 

c 

where N^. = rate at barium atoms are 

being emitted from the source, z is parallel -to B,and y 
is perpendicular to x and This expression is the 
cylindrical equivalent of the expression given in Eq. (88) 
for the spherical case. With the expression (96) for n, u , 

Eq. (87) can be integrated along the magnetic field to 
get the rate of change of the transverse distribu* ion in 
y of the field-line ion content. 


( 97 ) 


3N c (y,t) 

3t 


vNt 2 9 * / L ^ 

A - exp{-y^/R:(t)]e 1 


/? V(R c )R c (t) 


In this equation, the length of time, t*, that the barium 

atoms have Deen expos* d to ionizing radiation is not equal 

to t - t„ as * n the point-release case, but is equal to 

the time that it has taken the cloud to grow from radius 

R. to radius R„, 
c 

t*(t) = / 7^— — , (98) 

J R R(R’) 

i 

wh .*e R(R) = (dR/dx’ )V is gi' r en hi the venting model 
results discussed in Sec. IV. 


The uniform-in-x transverse profile of the field- 
line ion content, N (y), is obtained by integrating 
Eq. (97). It is convenient to change variables and integrate 
with respect to t* by making use of the relations 


dt * 

dt 


R(R C ) 


dt 


1 dx' 
V dR 


dR V 

C _ a 

dt ’ " V(R C ) 


(99) 


where the 
Eq. (98), 
With this 


first, second, and third equalities follow from 
the definition of R(R), and Eq. (95), respectively, 
substitution, ve obtain 
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( 100 ) 


N c <y) * 


™T f 

^ \ l 


expI-y 2 /R' 2 <t*) - t*/T 1 ) dt . 


R’(t*) 


where R'(t) * R[x*(t)] and x’ 


(t) ./' 


Vdt ’ is the axial 


position of a slab of released gas at time t after R * R i 

If we define a' as the radius of the wake at one 
time period after being exposed to ionizing radiation. 


a' - R’(x i ) . 


we may write Eq. (100) as 


fN„ 


N„(y) 


V v n a ' 

a 


i ' (y) 


v/here I'(y) is s. normalized field-line content; 

’exp - (y 2 /R' 2 (t ’ ' - t’/T. ] 


I ’ ( y ) * a’ 


/ 


R’(f ) 

The normalization is chosen so that 

I ’ ( y ) dy * /¥ a ’ 


dr 


0 i 


S' 


which is the same value that would result if I’(y)were 
exp(-y 2 /a’ 2 ) . 

In order to determine the properties of such 
a release we need to carry out the integral in Eq. (103). 
The venting-mode model determines the function R (t ) 




( 101 ) 


( 102 ) 


(103) 
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* 

once a specification of M q (which determines the charac- 
teristic length a) and R^ is made. For the Spruce ion 

7 

cloud, the detailed modeling resulted in a value for 

« 

the radius of the neutral cloud of 1.59 km = 1.26 a 

o 

when photoionization began. For this venting mode model, 

we take R. * 1.59 km. For the spherical point release 

f 1 

with a flattop distribution, a cloud radius equal to 

3 

1.26 a Q implies that the barium atom density, M q /(4ttR^/3) 

= p /2, was one-half the ambient atmospheric density. 

• a 

Vie apply this same criterion in the venting-mode model 
case , 

• M /VttR 2 -p/2, 

o' l a' 

which reduces to 

# VR 2 = 2V a 2 . 

l a 

Examination of the numerical integration of the 
venting mode equations indicates that criterion (104) 
for the nominal case occurs at R^ = 1.69 a. Thus, for 
the venting-mode equivalent of a Spruce-like point 
release at 190 km altitude with R.^ = 1.59 km, we have 
a = 0.94 km which corresponds to M = 9.7 kg/s. For 
these same conditions we find that a 1 = R’(t^) = 10.9 km. 

We have carried out the integral in Eq . (103) 
using these values, and have found that I'(y) is a bell- 


(104) 
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shaped curve that has a peak value I'(0) - 3.56 and has 
fallen to 0.57 = (l/e)I'(0) at y * 0.57 a’. Hence, a 
venting-mode release that produces barium atoms at the 
rate of 9.7 kg/s at 190 km altitude is calculated to 
produce a cylindrical trail of barium ions having a 
characteristic scale of a A = 0.57 a' = 6.2 km in the 
transverse direction and a maximum field-line content of 
1.56 flL/tV a')= 1.84 x 10 17 m -2 . 

1 a 

C. Scale size and magnitude of conductivity 

enhancements from '.kets and Spacelab 

In Secs. II, III, and IV we have discussed 
models that describe the expansion of neutral gas clouds 
released in an intermediate altitude regime in both a 
point release mode at both negligible and orbital 
velocity and f v--;n ting-mode at orbital velocity. In 
Secs. V-A and V-B we have described the application of 
the results of the neutral cloud modeling to the calculation 
of the deposition of ionization on magnetic field lines 
resulting from: a) a point 48-kg barium thermite 

release (M f = 3.56 kg) from a rocket at 190 km altitude; 
and b) a venting-mode barium thermite release = 9.7 kg/s) 

at orbital velocity at 190 km altitude. The procedures 
described for the former case produce results consistent 
with various detailed measurements made of such a release as 


06 


described in Ref. 7. The procedures described for the 
venting-mode case are derived from criteria establishe' 
in the point-release case, but there are no data available 
for comparison. 

In this section we derive the field-line-integrated 
Pedersen (FLIP) conductivity that results from the ionization 
enhancements that are created and describe how we model 
the characteristics of the conductivity modifications 
produced by releases of different magnitudes at different 
altitudes. 

The FLIP conductivity enhancement due to the 
deposition of ionization from barium thermite releases 
is defined by 


V h> • I 


/ 


e(z)n i (z) 
1 + e^(z) 


dz 


(105) 


where z is the coordinate along $, n^(z) is the distribution 
of ion concentration along the magnetic field line that 
passes through the peak value of n^z), e(z) = [Qt(z)]~* 
where fl * eB/m Ba is the barium ion gyrofrequency and t(z) 
is the altitude-dependent ion-neutral collision time, 
and h is the altitude of the peak in the ion concentration. 
Henceforth we will adopt a nominal value of ft « 35 s -1 
appropriate for a 5 * 10 - ^ T magnetic field. 



The ion-neutral collision time can be obtained 


from measurements of the mobility of barium ions in 
nitrogen gas** resulting in 

_ _ 8 . 6 x 10 15 

1 S 

n n 

where n is an effective neutral concentration in 
n 

3 

particles/m , 


n = 
n 


% 


* "o„ - °- 8 n c 


In deriving Eq. (106) it has been assumed that the collision 
cross section varies inversely as the square root of 
the temperature and that the collision cross section for 
oxygen molecules and atoms is the same as it is for 
nitrogen molecules. The fact that the lighter oxygen 
atoms are less efficient in stopping a heavier barium 
ion is expressed by the appropriate Langevin factor, 

O + m Ba/ m N ^* + m Ba^ m 0^ * ~ 0,8 

A 


where m XT 
N 2 

and oxygen 


and mg are the masses of the nitrogen molecules 
atom respectively. 


The distribution of ionization along the magnetic 
field through the peak value is governed by Eq. (87) with 
the right-hand-side evaluated at r = 0. As discussed in 


(106) 


(107) 
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I 


9 


» 


f 


I 


9 


1 


Sec. V-A, after the passage of several time periods 
(typically of order 1 minute after release) most of the 
ionization has been deposited on field lines. The 
extent of the distribution of ionization along field 
lines is still small at that time compared to the atmos- 
pheric neutral-density scale-height H . Hence, the value 

n 

of the FLIP conductivity in Eq. (105) resulting from the 
deposition of ionization on field lines can be approximated 
by 


Vv * I 


1 + e 


2 N c ~ 3 ‘ 2 


X 10 


-15 


N 


mho 


1 + e‘ 


(108) 


where e r is evaluated at the altitude of release, h r , id 

-2 

N c is the peak field-line content in m through the 
cen’.er of the cloud given by 


N c = N c (0.») 


fN T 

1.65 

7t a 


(109) 


for point releases where N c (r,t) is defined by Eq . (91), and 

fN T 

N - N (0) = 1.56 — — (110) 

c c V /if a' 

a 

for venting-mode releases where N p (y) is deiined by Eq . (102). 

Note that N^,/V a is the number of barium atoms that the 
source emits per unit length and f is the fraction of 
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these atoms that are phctoionized. 

The numerical coefficients in Eqs. (109) 
and (110) were derived on the basis of calculations 
appropriate for the nominal 190 km altitude releases 
described earlier (M q * 3.56 kg and M q = 9.7 kg/s, 
respectively). We have not carried out the detailed 
ion deposition calculations indicated by Eqs. (92) 
and (103) for releases of different amounts of barium 
vapor at different altitudes. Instead, we approximate 
the field-line content produced by other releases by 
the expressions (109) and (110). In these expressions, 

N t and N,p depend directly on the amount of released 
vapor or the rate of vapor release, respectively, and 
f, given by Eq. (85), depends on the altitude of 
release. The scale sizes of the neutral clouds one 
time period after photoionization, a and a', depend 
in a complicated way on both factors as described 
below. 

By combining Eqs. (109) and (110) with Eq. (108), 
we obtain expressions for the FLIP conductivity as a 
function of amount of released gas and altitude, 

_ 91 e r fN T 

2 (u r ) - 1.68 x lO - ^ 1 mho , 


(in) 
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mho 


( 112 ) 



U r> 


2.82 x 10 -21 


fN„ 


1 V a a' 


for the cases of point releases, and venting-mode releases, 
respectively, with a, a', and V expressed in km, km, and 

c* 

km/s , respectively. N_ and are the number of barium atoms 
and number of barium atoms per second, respectively. 


Equation (27) in Sec. II-B provides an excellent 
analytic approximation to the neutral cloud radius, R(t) 
for R > a Q . For the cas- of stationary point releases, 
we take for the characteristic size of the neutral cloud 
at one time period, t., after photoionization commences, 
the expression 


a * 1.2 a 

o 


Z 1 - 15 T i 

\ a o 



0.8 


~\ 1/2 


0.44 


km 


(113) 


where a Q is in km, is given by Eq. (86), and S i is a 
parameter reflecting the time delay before the onset of 
photoionization. V»e assume that photoionization commences 
when the optical depth of the cloud to the metastable- 
exciting radiation approaches one. This assumption is 
equivalent to a critical value for the product of the 
barium density times the cloud radius or 
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(114) 


P i R i 


3M 

o 

4ir R 2 


0.34 kg /km 2 


where the numerical value was obtained from the Spruce 

case. If we use Eq. (27) with S. = ft t./a^ in oraer to 

1 O 1 o 

evaluate Eq. (114), we obtain 


S. = (0.44 + 2p a J 1 ' 25 (115) 

1 a O 

o 

where pa is expressed in kg /km . This value used in 
a o 

Eq. (113) determines the dependence of a on mass of 
released gas and altitude. 


we take 


a 


For the case of a point release at orbital velocity, 


2.68 a 

o 


f 1 - 15 T i 

\ a o 



T 1/2 


0.44 


km 


(116) 


which is larger than Eq. (113) by about a factor of 2.3 . 

Ag .in by applying the criterion (114), we obtain 

S. * (0.44 + 0.4 pa) 1-25 . (117) 

1 a O 

The expressions (115) and (117) are not accurate whenever 
the corresponding critical radius R^ < a Q . This condition 
occurs only atove a critical altitude depending on M q when 
" 0.6 and is much smaller than the term 1.15 t /a to 
which it is added. Hence this approximation to a does not 
lead to significant error. 
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In the case of a venting-mode release, Eq (76) 


can be integrated in order to obtain time as a function 
of wake radius which provides an implicit R(t): 

K 2 V (t - t ) = R(R 2 /a 2 + B)* + a(2 - B) in [R/a + (R 2 /a 2 + B)*] 

2- O 

where t is an arbitrary constant of integration. This 
equation cannot be inverted to obtain R(t) implicitly but 
we have found that the expiession 


R( t ) = a 


K 2 V 


a 


. . 0 


t + S i + 0.2 - in 


11/2 


with a in km and K V = 12.2 km/s provides an adequate 

representation for the radius for the nominal case for 

R > 1.2 a with = 1.24. Because d.nd t are additive, 

the time t in this equation can be interpreted as the time 

lapse after photoionization commences i.f we scale 3^ 

properly with pa in order to reflect the critical time 
a 

to reach a critical radius R^ . The criterion (ll- 1 ) in 
this case reduces to 


U 


p . R. 

l l 


V R i 


° =0.34 kg/km 2 


or 


V.R i R 2 /a 

V a a (R 2 + 2.8 a 2 )^ 


oP a a 


(113) 


(119) 


(120) 


( 121 ) 
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) 


o 

with p a expressed in kg/km . The middle expression is 

cL 

an approximate expression which is reasonably accurate 

for R. > 1.2 a. By solving this quadratic equation for 
2 

(R./a) in terms of pa, using this value in Eq. (118) 

1 a 

and adjusting t so that S. =4.05 when R./a * 1.69 as 
for Spruce, we obtain the following expression for 
for the venting mode case, 



Equation (119) then defines the characteristic size of 
the neutral cloud wake as a' = £( 1 ^. 

As discussed in Sec. V-A, the transverse scale 
size of the ion cloud field-line content defined as the 
half-width a (1/e) of the maximum value is given as 

a A = 0.7 a 

in the case of point, releases and 

a x =0.57 

in the case of venting-mode releases. 

Figures 21-26 show contours of constant value; 
of che transverse seal' size, a^, and FLIP conductivity, 
£ p , for the three release mode cases. The shaded portions 


r* 


( 122 ) 


(123) 


024) 


94 


of these figures represent regions in which the snowplow 
models used, applicable only to the intermediate altitude 
regime, are not valid. The contours should not be used 
in the shaded regions because their actual form is 
qualitatively different from those shown. Fortunately, 
the allowed regions correspond to regions of many 
interesting suggested experiments. 

As a typical result, a 48 kg barium release at 
190 km altitude will produce a FLIP conductivity of approx- 
imately 30 mho, several times the ambient ionospheric 
conductivity, and will have a transverse scale size just 
less than 3 km. The extra energy available in the shuttle 
release causes the gas to expand faster resulting in the 
barium ionization being spread over a larger area. The 
correspondingly lower density of ionization results in 
lower conductivity. If the same amount of chemical 
payload is released from shuttle, the resulting conductivity 
is lover by a factor of 5 from that produced by a release 
from a rocket and the transverse size of the ion cloud 
is increased by a factor of approximately 2.3. 

As time passes, th* conditions envisioned by 
tuis modeling change and the distribution of ionization 
both parallel and perpendicular to § chant . . After 
about a minute following release the ionization clouds 
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that are produced have the characteristics that we have 
described. Beyond this time, the distribution c f ionization 
along the magnetic field is governed by Eq. (87) with the 
right-hand side equal to zero. The z-component of the 
ionization velocity is given by 


1Z 


= V_._ + g„T +k(.T e + T i )T — 
w ' z 


nz 




(125) 


where v and g„ are the components of neutral wi*... 
nz z 

velocity and gravity parallel to 2, respectively, and 
T^ are the electron and ion temperatures, respectively, 
and x is the ion-neutral collision time given by Eq. (106). 

The resulting evolution of the distx'ibut ion parallel to ^ 
as a result of neutral wind, gravity and pressure gradients 
in an exponentially varying atmosphere whose effective 
number density, n, defined in Eq, (107), varies as 

n(h) = n(h r ) exp[-(h - h r )/H n i (126) 

has been discussed i' detail in Sec. 4 of Ref. 7. 

Initially the distribution is Gaussian, then the 

peak density descends in altitude and the distribution 

in altitude evolves into a Chapman-like ,.ayer with the 

density falling off rapidly on the underside and 

exponentially on the top side. An analytical solution o + ' 

Eq. (87; for n.(z,t) with v. given by Eq . (1*25) ha? been obtained 
1 1 z 


ft*? 
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f 


that describes the evolution as a function of time for 
various models of altitude-dependent neutral winds. By 
analytically integrating the solution, the remarkable 
rer ult was obtained that at a later time when the peak 
in rhe ionization density has descended to altitude h, 
the FLIP conductivity has the simple value 

2 p (h) ” 2 p (h r )exp[-(h - h r )/H n ) . 

We remark that this fully-time-dependent result obtained 
in Ref. 7 does not require that z be in the vertical 
direction, and that the winds can be either constant or 
altitude dependent with an exponential profile. 

In the direction transverse to the magnetic 

field, the situation is far more complex. A neutral 

wind blowing across a plasma gradient leads to the 

12 

well-known gradient-drift instability. The result is 
that the coherent picture of a conductivity modification 
no longer applies because the ionization striates. The 
typical time scale for point releases from rockets between 
150 and 250 km for striations to onset is 5 to 20 minutes. 
Simple theoretical modeling indicates that this onset 
time scales with £a /v where ? is a factor that depends 
on the magnitude of the conductivity enhancement and 
typically lies between 5 and 20. Hence, concepts to 
affect ionospheric and magnetospheric current flows by 

97 


( 127 ) 



conductivity enhancements such as we have been describing 
can probably only b» planned for the first 10 minutes 
following release. 
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Fig re Captions 


FIG. 1. Cloud radius vs. time for various values 

of the ratio of specific heats, y, and density profile 
parameter, a. Translational velocity at release is 7.8 
km/s except for the curve u"0. 

FIG. 2. Temperature, radius, and distance traveled 

vs. time for a 100 kg gas release at orbital velocity at 
200 km altitude. 

FIG. 3. The evolution of the normalized density 

distribution, g(z, t d ), from an initial flat-top to a 
Gaussian as a result of diffusion for several values of 
the diffusion parameter ip. 

FIG. 4. The evolution of the normalized density 

distribution, g(z, t^), from an initial shell to a 
Gaussian as a result of diffusion for several values of 
the diffusion parameter t^. 

FIG. 5. Normalized diffusion function 7 (s) for a 

gas with y«5/3 in terms of which the diffusion parameter 
ip is defined according to Eq. (38). The dashed curves 
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show the values of t(s) if the snowplow were stopped at 
s«18 and 46 respectively. 

FIG. 6. Number density, n, mass density, p ft ,and 

mean-f ree-path , A , as a function of altitude based on the 
CIRA 1972 model atmosphere. Note the logarithmic altitude 
scale. The solid curves correspond to an exospheric 
temperature * 1200°; the dashed curves correspond to 

©X 

T - 700° and 2000°. 

x? A 

FIG. 7. Contours of constant characteristic length, 

a, for a CIRA 1972 model atmosphere with T ex » 1200° in 
the altitude versus mass-venting-rate plane. The dashed 
curves show the shape of the contours for T ■ 700° and 
2000 °. 


FIG. 8. Contours of constant a/ A in the altitude 

versus mass-venting-rate plane. The dashed curves show 
the shape of the contours for T_„ ■ 700° and 2000°. 

©X 

FIG. 9. Radial expansion velocity as a function 

of normalized radius for different density profiles and 
specific internal energies. 



FIG. 10. K&di&l expansion velocity as a function 

of normalized radius for different assumed initial 
velocities. 

FIG. 11. Radial expansion velocity as a function 

of normalized radius for different values of the ratio 
of specific heats. 

FIG. 12. Temperature, defined according to Eq. (56), 

as a function of normalized radius for different assumed 
initial expansion velocities and ambient atmospheric 
temperature. 

FIG. 13. Temperature, defined according to Eq. (56), 

as a function of normalized radius for different values 
of the ratio of specific heats. 

FIG. 14. Axial velocity in the rest frame of the 

canister as a function of normalized radius for different 
assumed initial expansion velocities. 

FIG. 15. Normalized radius versus normalized axial 

position for various values of gas release parameters 
plotted on a logarithmic scale. 
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FIG. 16. Normalized radius versus normalized axial 

position for various values of gas release parameters plotted 
on a linear scale. 

FIG. 17. Normalized radius versus normalized axial 

position for different assumed initial expansion velocities 
plotted on a logarithmic scale. 

FIG. 18. Normalized travel distance versus normalized 

axial position for different density profiles and specific 
internal energies. 

FIG. 19. Normalized travel distance versus normalized 

axial position for different assumed initial expansion 
velocities. 

FIG. 20. Normalized travel distance versus normalized 

axial position for different values of the ratio of 
specific heats. 

FIG. 21. Contours of constant transverse scale size 

of the conductivity modification in the altitude versus 
chemical payload weight plane for point barium releases 
from rockets. In this and the next 5 figures the model is 
invalid in the shaded regions. 
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FIG. 22. Contours of constant f ield-lice-integrated 

Pedersen conductivity in the altitude versus chemical 
payload weight plane for point barium releases from 
rockets . 

FIG. 23. Contours of constant transverse scale 

size of the conductivity modification in the altitude 
versus chemical payload weight plane for point barium 
releases at orbital velocity. 

FIG. 24. Contours of constant field-line-integrated 

Pedersen conductivity in the altitude versus chemical payload 
weight plane for point barium releases at orbital velocity. 

FIG. 25. Contours of constant transverse scale size 

of the conductivity modification in the altitude versus 
chemical payload release-rate plane for venting mode 
releases at orbital velocitj». 

FIG. 26. Contours of constant field-line-integrated 

Pedersen conductivity in the altitude versus chemical 
payload release-rate plane for venting mode releases 
at orbital velocity. 
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FIG, l. Cloud radius vs. time for various values of the ratio 
of specific heats, y, and density profile parameter , a. 
Translational velocity at release is 7.8 km/s except for the 


TEMPERATURE. RADIUS, DISTANCE TRAVELED VERSUS TIME 



FIG. 2. Temperature, radius, and distance traveled vs. time 
for a 100 ktf pas release at orbital velocity at 200 km altitude 
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FIG. 3. The evolution f the normalized density distribution, g(z, 
from an initial flat-top to a Gaussian as a result of diffusion for 
several values of the diffusion parameter T n . 
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FIG. 6 . Number density, n. mass density, p , and 
mean- f ree-path , X, as a function of altiludU based 
on the CIRA 1972 model atmosphere. Note the lop" - 
rithmic altitude scale. The solid curves corres- 
pond to an exospheric temperature T ■ 1200°; the 
dashed curves correspond to ■ 70$° and 2000°. 
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FIG. 7. Contours of constant characteristic length, 
a. for a CIRA 1972 model atmosphere with T » 1200° 
in the altitude versus mass-vent ing-rate pfane. The 
dashed curves show the shape of the contours for T 
■ 700° and 2000°. ex 
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FIG. 9. Radial expansion velocity as a function of 
normalized radius for different density profiles and 
specific internal energies. 
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FIG. 10. Radial expansion velocity as a function of 
normalized radius for different assumed initial velocities, 
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FIG. 11. Radial expansion velocity as a function of 
normalized radius for different values of the ratio o 
specific heats. 
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FIG. 12. Temperature, defined according to Eq. (56), as 
a function of normalized radius for different assumed 
initial expansion velocities and ambient atmospheric 
temperature . 
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FIG. 13. Temperature, defined according to Eq. (56), as 
a function of normalized radius for different values of 
the ratio of specific heats. 
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FIG. 14. Axial velocity in the rest frame of the canister 
as a function of normalized radius for different assumed 
initial expansion velocities. 




119 


I 


NOMINAL CASE; y « 1.4, a - 0.5. 

C n * 0.89 KM/S 


R 0 -A 
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FIG. 15. Normalized radius versus normalized axial 
position for various values of gas release parameters 
plotted on a logarithmic scale. 
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FIG. 16. Normalized radius versus normalized axial 
position for various values of gas release parameters 
plotted on a linear scale. 
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FIG. 17. Normalized radius versus normalized axial 
position for different assumed initial expansion 
velocities plotted on a logarithmic scale. 


122 



L/A 





9 

t 


X/A 

FIG. 18. Normalized travel distance versus normalized 
axial position for different density profiles and specific 
internal eneigies. 
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FIG, 19. Normalized travel distance versus normalized 
axial position for different assumed initial expansion 
velocit ies . 
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FIG. 20. Norm.'. H zed travel distance versus normalised 
axial position for different values of the ratio of 
specific heats. 
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FIG. 21. Contours of constant transverse scale 
size of the conductivity modification in the 
altitude versus chemical payload weight plane for 
point barium releases from rockets. In this and 
the next 5 figures the model is invalid in the 
shaded regions. 
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FIG. 22. Contours of constant field-line-integrated 
Pedersen conductivity in the altitude versus chemical 
payload weight plane for point, barium releases from 
rockets . 
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FIG. 2A. Contours of constant field-1 ine-integrated 
Pedersen conductivity in the altitude versus 
chemical payload weight plane for point barium 
releases at orbital velocity. 
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FIG. 25. Contours of constant transverse scale 
size of the conductivity modification in *he 
altitude versus chemical payload release-rate 
plane for venting mode releases at orbital 
velocity. 
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FIG. 26. Contours of constant field-line- 
integrated Pedersen conductivity in the altitude 
versus chemical payload release-rate plane for 
t venting mode releases at orbital velocity. 
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